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EUSTIS  DIRECTORATE  POSITION  STATEMENT 


A method  for  determining  rotor  hub  load*  (force*  and  moment*)  directly  from  measured  fuselage  accelerations  was 
successfully  demonstrated.  The  method,  called  Force  Determination,  was  verified  on  a laboratory  modal  repre- 
sentative of  a helicopter  in  term*  of  mats  and  stiffness  distributions,  natural  frequency,  and  structural  damping. 
Following  calibration  of  the  5-foot  dynamic  model,  known  vertical  and  lateral  hub  forces  were  applied  simulta- 
neously and  the  fuselage  responses  recorded.  Using  these  responses  and  the  calibration  data,  the  hub  forces  were 
calculated  using  Force  Determination.  The  overall  accuracy  was  excellent.  In  general,  the  errors  in  the  predicted 
hub  forces  were  within  the  achievable  accuracy  limits  of  the  instrumentation.  The  errors  were  on  the  order  of  2-5% 
in  amplitude  and  2-4  degrees  in  phasing.  In  cases  where  the  model  was  calibrated  at  a nominal  gross  weight  and  eg 
and  then  tested  with  force  predictions  made  at  a condition  of  extreme  grots  weight  or  eg,  errort  of  22-24%  resulted. 
However,  when  calibrated  for  the  tested  gross  weight  and/or  eg,  the  errort  returned  to  the  above-cited  excellent 
results.  This  suggests  that  for  actual  aircraft  use,  calibration  of  the  vehicle  should  not  simply  be  done  once  but  for 
a set  of  conditions  representing  the  operational  gross  weight  and  eg  range. 

It  has  been  concluded  that  Force  Determination  is  a technique  "in-the-bud"  which,  if  successfully  developed  and 
implemented,  could  "blossom"  into  a technological  breakthrough  with  far-reaching  developmental  and  operational 
possibilities.  The  possibilities  foreseen  are:  (1)  Force  Determination  could  expedite  helicopter  development. 

During  flight  test  development  hub  loads  over  the  mission  profile  could  be  determined.  In  turn,  engineers  could 
duplicate  these  loads  on  an  actual  helicopter  in  an  accelerated,  ground-based,  simulated  flight  test  to  establish 
system  and  component  reliability.  Failures  and  their  modes  would  be  exposed  much  earlier  in  the  development 
phase.  Being  ground-based,  conditions  pertinent  to  damage/failure  modes  could  be  easily  observed,  permitting 
timely,  effective  redesign  prior  to  fleet  introduction.  And  (2),  operationally,  "lead-the-fleet"  endurance  testing 
is  foreseen.  This  would  allow  low-cost,  highly-responsive  reliability  testing  for  detecting  problems  before  they  arise 
in  the  field  and  allow  the  evaluation  of  corrective  redesigns  before  they  are  deployed  in  the  -field.  With  "lead-the- 
fleet”  testing,  it  it  envisioned  that  sensors,  small  units  not  encumbering  vehicle  utility,  would  be  "planted"  in  air- 
craft delivered  to  the  fleet  to  record  in-flight  responses;  and  hub  loads,  determined  from  these  responses,  could  be 
used  for  accelerated  endurance  testing  to  induce  failures  before  they  occur  in  the  fleet. 

Demonstration  of  this  new  method  on  an  Army  inventory  helicopter  is  planned. 

This  program  was  conducted  under  the  technical  cognizance  of  Joseph  H.  McGarvey,  Military  Operations 
Technology  Division. 
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Results  of  laboratory  tests  on  a helicopter  dynamic  model  show  that  the  external  vibratory 
forces,  such  as  those  of  main  and  tail  rotors,  acting  on  a helicopter  can  be  accurately  de- 
termined from  flight  tests  using  only  accelerometers  in  the  fuselage,  r 


The  contract  was  in  three  tasks:  (I)  Design  and  fabrication  of  a model  representative  of 
helicopters  in  terms  of  mass  and  stiffness  distribution  as  well  as  natural  frequency  and 
structural  damping;  (II)  Dynamically  calibrate  the  model  followed  by  exciting  the  hub  with 
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known  forces  and  verify  the  method  of  Hub  Force  Determination  by  calculating  what  the 
forces  must  have  been,  using  the  measured  fuselage  accelerations  and  the  dynamic  calibra- 
tion data  as  a basis;  and  (Hi)  Dynamically  calibrate  the  model  by  both  direct  hub  exci- 
tation and  single-point  excitation  at  any  other  convenient  point.  This  *hird  task  was  to 
demonstrate,  that  single-poinf'excltation  can  yield  fuselage  accelerations  which  together  with 
the  methods  of  Modal  Acceleration  testing  and  force  Determination  could  be  used  as  a 
viable  structural  dynamic  test  method,  yielding  the  same  results  as  if  the  structure  had  been 
excited  at  the  hub.  If  successful,  a considerable  time  and  cost  savings  would  be  realized. 

It  would  no  longer  be  necessary  to  shake  the  helicopter  at  the  hub.  Instead,  shaking  could 
take  place  with  the  excitation  applied  at  any  convenient  point  other  than  a node. 


Following  calibration  of  the  5-foot  dynamic  model,  known  vertical  and  lateral  hub  farces 
were  applied  simultaneously.  Different  flight  speeds  and  maneuvers  were  represented Y>y 
different  magnitudes  and  phasings  of  the  shaker  forces  at  the  hub.  This  was  done  for 
various  "rotor  speeds"  (shaker  frequencies)  as  well  as  a broad  range  of  gross  weight  and 
CGs.  For  each  flight  condition  the  amplitudes  and  phasings  of  the  twelve  fuselage 
accelerometers  were  recorded.  Using  these  recorded  responses  and  the  calibration  data, 
the  hub  forces  were  calculated  using  Force  Determination.  A total  of  81  test  cases  were 
run  covering  54  combinations  of  vertical  and  lateral  forces  over  phase  angles  ranging  from 
0°  to  180°.  The  overall  accuracy  of  Force  Determination  was  shown  to  be  excellent.  In 
general,  the  errors  in  the  predicted  hub  forces  were  within  the  achievable  accuracy  limits 
of  the  instrumentation.  The  errors  were  on  the  order  of  2-5%  in  amplitude  and  2-4 
degrees  in  phasing.  In  cases  where  the  model  was  calibrated  at  a nominal  gross  weight 
and  eg  and  then  tested  with  force  predictions  made  at  a condition  of  extreme  gross  weight 
or  eg,  errors  of  22-24%  resulted.  However,  when  calibrated  for  the  tested  gross  weight 
and/or  eg,  the  errors  diminished  to  the  above-cited  excellent  results.  This  suggests  that 
for  actual  aircraft  use,  calibration  of  the  vehicle  should  not  simply  be  done  once  but  for 
a set  of  conditions  representing  the  operational  gross  weight  and  eg  range.  These  results 
were  equally  good  whether  the  model  was  actually  shaken  at  the  hub  or  at  some  convenient 
fuselage  station.  Thus,  the  additional  advantages  of  single-point  excitation  appear  to  be 
realizable. 
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Langley  Directorate,  USAAMRDL. 

Force  determination  on  the  Bell  UH-1H  DAVI-equipped  aircraft  was  done 
under  the  technical  supervision  of  Robert  Jones.  Alex  Berman  served 
as  analytical  consultant  and  Nicholas  Giansante  as  experimental  con- 
sultant in  the  investigations.  The  helicopter  dynamic  model  was 
constructed  by  Harry  A.  Cooke  and  W.  T.  Rhodes.  Among  those  contributing 
to  the  instrumentation  work  were  D.  Stewart,  I.  N.  Johnson,  and 
R.  F.  Metzger.  G.  T.  Neu  was  engineering  editor  of  the  final  report. 
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Reliability  and  Vibration 

The  Army  has  a continuing  goal  of  improving  the  reliability  of  its  air- 
craft and  particularly  its  helicopters.  Although  in  recent  years 
availability  rates  of  helicopters  have  approached  and  sometimes  exceeded 
those  of  fixed-wing  aircraft,  it  has  only  been  at  a very  high  cost  in 
maintenance  man-hours,  spares  bought  and  stocked,  and  excessive  depot 
overhaul  of  aircraft  and  components. 

A fundamental  cause  of  the  high  maintenance  man-hours,  the  frequent 
component  replacement,  and  the  generally  low  helicopter  reliability, 
is  the  high  vibratory  loads  the  machine  experiences.  In  general,  the 
highest  vibratory  loads  are  generated  by  the  rotor  and  occur  at  blade 
passage  frequencies,  although  there  are  significant  loads  at  other 
frequencies.  All  major  helicopter  manufacturers  have  published  technical 
reports  and  papers  outlining  their  efforts  to  reduce  these  vibrations 
and  hence  increase  reliability.  One  study  (Reference  1),  performed  by 
Sikorsky  under  contract  to  the  Eustis  Directorate  of  USAAMRDL,  found  very 
large  increases  in  helicopter  reliability  resulting  from  the  reduction 
in  blade  passage  frequency  forces  transmitted  to  the  fuselage. 

There  are  three  major  technical  thrusts  in  the  industry  for  increasing 
the  reliability  of  helicopters  by  reducing  the  vibratory  forces  and 
moments  transmitted  to  the  fuselage  from  the  rotor:  (1)  improved  rotor 
design;  (2)  rotating  system  dynamic  absorbers  (hub  absorbers,  blade 
absorbers);  and  (3)  rotor  isolation  (anti resonant,  conventional,  and 
active).  Excellent  progress  has  been  made  in  these  areas  in  recent 
years  in  spite  of  technical  problems. 

Force  Determination 

A technique  developed  by  Kama n Aerospace  Corporation  under  contract  to 
the  Eustis  Directorate,  USAAMRDL,  is  a solution  to  many  of  the  technical 
problems  which  frustrate  continued  progress  in  the  three  areas  listed 
above.2  This  technique  is  called  force  determination.  Force  determina- 
tion  is  a method  for  determining  the  magnitudes  and  phasings  of  vibratory 

’Angelo  C.  Veca,  VIBRATION  EFFECTS  ON  HELICOPTER  RELIABILITY  AND  MAINTAIN- 
ABILITY, Sikorsky  Division,  United  Technologies  Corp. ; USAAMRDL  Technical 
Report  73-11,  Eustis  Directorate,  U.  S.  Army  Air  Mobility  Research  and 
Development  Laboratory,  Fort  Eustis,  VA,  April  1973,  AD  766307. 

2W.G.  Flannel ly  and  N.  Giansante,  EXPERIMENTAL  VERIFICATION  OF  SYSTEM 
IDENTIFICATION,  Kaman  Aerospace  Corporation;  USAAMRDL  Technical  Report 
74-64,  Eustis  Directorate,  U.  S.  Army  Air  Mobility  Research  and  Develop- 
ment Laboratory,  Fort  Eustis,  VA,  August  1974,  AD  784190. 
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rotor  forces  and  moments  (and  other  external  forces  and  moments)  on  a 
helicopter  in  flight  through  accelerometer  measurements  on  the  fuselage 
only. 

Rotor  Analysis  Correlation  and  Vibration  Effects  of  Proposed  Rotors 

In  1974,  Dr.  R.  A.  Ormiston  of  USAAMRDL  conducted  a comparative  test  of 
rotor  loads  analysis  methods  of  the  U.  S.  military  helicopter  industry.3 
"Large"  to  "very  large"  differences  in  root  shears  were  found  among  the 
analytical  methods  of  the  industry  but  experimental  data  were  not  avail- 
able for  establishing  the  accuracy  of  any  loads  prediction  method. 

Dr.  Ormiston's  findings  raise  disquieting  questions  on  the  validity  of 
rotor  model  shear  predictions. 

A great  deal  of  effort  in  the  helicopter  industry  is  devoted  to  improve- 
ments in  the  mathematical  analysis  of  rotors.  If  analysts  had  at  hand 
accurate  in-flight  measurements  of  the  magnitudes  and  phasings  of  hub 
shears  and  moments  in  various  maneuvers,  they  could  then  cheaply  and 
conveniently  correlate  predictions  from  their  mathematical  rotor  models 
over  a range  of  forward  speeds.  This  would  enable  them  to  select  more 
rapidly  the  most  suitable  mathematical  models  for  each  case  and  thereby 
permit  faster  and  more  effective  improvements  in  rotor  analysis  computer 
programs. 

Although  prediction  of  hub  shears  and  moments  is  not  the  sole  measure  of 
rotor  model  accuracy,  force  determination  allows  these  predictions  to  be 
checked  against  known  flight  loads.  This  check  serves  as  a valuable 
indication  of  the  model's  overall  validity  and  allows  the  engineer  to 
weed  out  those  models  demonstrably  unsuitable  or  select  from  among 
alternative  models. 

Having  a rotor  model  whose  validity  has  been  verified  within  the  flight 
regimes  of  interest  will  be  of  inestimable  benefit.  New  rotors  proposed 
for  existing  aircraft  can  be  analyzed  to  evaluate  the  fatigue-producing 
vibratory  loads.  If  the  analysis  shows  improvement,  the  manager  can 
move  from  design  to  hardware  with  much  more  confidence  than  heretofore 
possible.  In  short,  force  determination  will  permit  the  verification  of 
mathematical  formulations  of  rotors  based  on  information  gathered  in 
flight  in  a manner  that  has  not  been  possible  previously. 

Rotor  Absorber  and  Isolator  Design 

The  motions  which  a rotor  absorber  or  a rotor  isolator  must  undergo,  and 
therefore  the  size  of  such  devices,  depend  on  the  magnitudes  and  phasings 
of  rotor  shears  and  moments  to  which  the  absorber  is  subjected,  and  on 


3 R.  A.  Ormiston,  COMPARISON  OF  SEVERAL  METHODS  FOR  PREDICTING  LOADS  ON 
A HYPOTHETICAL  HELICOPTER  ROTOR,  Proceedings  of  the  Meeting  of 
Specialists  on  Rotorcraft  Dynamics,  American  Helicopter  Society  and 
Ames  Research  Center,  Sponsors,  13-15  February  1974. 


the  mobility  of  the  fuselage  at  its  attachment  points.  The  fuselage 
mobility  also  determines  bandwidth  for  a given  weight.  Force  determina- 
tion provides  both  pieces  of  necessary  information  - excitation  and 
mobility  - from  actual  testing,  thereby  eliminating  much  of  the  trial  and 
error  in  flight  testing  rotor  isolator  and  absorber  designs. 

Rotor  absorbers  and  isolators  do  not  usually  reduce  all  the  rotor  forces 
and  moments  transmitted  to  the  fuselage,  and  even  improved  rotor  designs 
may  decrease  some  of  the  hub  loads  while  increasing  others.  The  engineer 
needs  to  know  the  phasing  among  the  hub  loads  to  decide  which  should  be 
reduced  to  provide  the  lowest  vibration  at  any  given  point  on  the  ship. 
Reducing  only  the  largest  load  may  actually  increase  the  total  vibration 
at  critical  fuselage  points,  as  V.  Sankewitsch  of  Boeing  Vertol  once 
pointed  out  with,  in  essence,  the  following  physical  explanation: 

The  resultant  vibration,  say  in  g's,  at  any  point  on  the  fuselage  is  the 
resultant  of  the  vector  sum  of  the  vibrations  at  that  point  due  to  each 
of  the  hub  excitations.  An  accelerometer  on  the  centerline  of  the  heli- 
copter shows  the  resultant  of  the  accelerations  due  primarily  to  vertical 
hub  shear,  longitudinal  hub  shear,  and  hub  pitching  moment  as  shown 
hypothetically  in  Figure  1.  If  the  longitudinal  force,  which  has  the 
largest  effect  in  this  hypothetical  example,  is  eliminated,  the  resultant 
vibration  gets  significantly  worse. 
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Removing  only  the  largest  hub  force  or  moment  can  actually  make  the 
vibration  worse  because  of  the  phasing  of  mobility  and  loads.  Force 
determination  testing  would  provide  all  the  necessary  decision-making 
information. 

Rotor  Mass  in  Shake  Testing 

With  the  results  of  force  determination,  flying,  and  the  dynamic  calibra- 
tion in  force  determination  ground  testing,  the  engineer  could  predict 
the  effects  of  changing  any  of  the  transmitted  hub  loads. 

Shake  tests  of  helicopters  have  had  the  unknown  rotor  impedances,  which 
change  the  natural  frequencies,  "simulated"  with  purely  conjectural 
lumped  masses  (sometimes  equal  to  the  rotor  mass  and  sometimes  half  the 
rotor  mass,  depending  on  the  preference  of  the  individual  engineer).  The 
erroneous  idea  that  some  lumped  mass  is  a practical  approximation  to  rotor 
impedance  probably  originated  with  a misinterpretation  of  Coleman's 
pioneering  work  of  the  1940's  in  mechanical  instability,  also  known  in- 
accurately but  popularly  as  "ground  resonance"  from  the  early  autogiro 
accidents.4  Coleman  used  a lumped  rotor  mass  at  the  hub  for  testing  as 
a convenience  because  he  approximated  hub  modal  mass  by  shifts  in  the 
natural  frequency  with  hub  mass  changes,  a technique  now  known  to  be 
accurate  only  in  a single  degree  of  freedom  system.  However,  Coleman 
accounted  for  the  lumped  mass  in  his  analysis  by  coupling  the  hub  response 
of  the  grounded  fuselage  with  his  in-plane  rotor  analysis  of  rigid  hinged 
blades.  He  did  not  simulate  hub  impedance  by  a lumped  mass. 

A study  of  the  work  of  Coleman  and  all  that  has  followed  raises  the 
serious  question  of  whether  the  idea  of  a total  helicopter  natural 
frequency  has  any  definable  meaning. 

In  recent  years.  Dr.  D.  R.  Gaukroger  and  others  of  the  Royal  Aeronautical 
Establishment  at  Farnborough,  England,  have  succeeded  in  making  actual 
measurements  of  the  impedances  of  a model  helicopter  rotor  for  the  first 
time. 5 Their  experiment's  results  clearly  show  that  the  spectral  re- 
sponse of  a rotor  is  not  at  all  approximated  by  a mass  or  other  simple 
impedance  element.  This  particular  conclusion  should  not  be  surprising 
as  it  has  been  known  since  1942  that  the  fixed  system  equations  of  even 
the  simple  in-plane  case  of  articulated  rigid  blades  have  periodic  co- 
efficients. The  existence  of  periodic  coefficients  makes  it  impossible 


R.P.  Coleman  and  A.M.  Feingold,  THEORY  OF  SELF-EXCITED  MECHANICAL 
OSCILLATIONS  OF  HELICOPTER  ROTORS  WITH  HINGED  BLADES,  National  Advisory 
Committee  for  Aeronautics  Technical  Note  3844,  February  1957. 

R.  Cransdale,  D.R.  Gaukroger,  and  C.W.  Skingle,  A TECHNIQUE  FOR 
MEASURING  IMPEDANCES  OF  A SPINNING  MODEL  ROTOR,  Royal  Aircraft 
Establishment  Technical  Report  71092;  Ministry  of  Defence,  Farnborough 
Hants,  England,  May  1971. 
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to  rigorously  define  a "lambda  matrix"  with  roots  (i.e.,  eigenvalues  or 
characteristic  values)  so  that,  according  to  the  Hamilton -Cay ley  Theorem 
(circa  1822),  "natural  frequencies"  can  be  mathematically  defined  in  the 
form  in  which  we  would  like  to  use  them  in  helicopter  engineering.  It 
seems  reasonable  to  expect  that  no  solutions  to  the  problem  will  be  found 
before  the  problem  is  defined. 

DEFINING  THE  PROBLEM  - FORCE  DETERMINATION 

What  an  engineer  needs  to  know  to  make  a helicopter  shake  test  is  not 
some  value  of  the  dubious  abstract  concept  of  "natural  frequency"  of  the 
total  system  but  the  forces  and  the  phasing  of  the  excitations,  how  the 
helicopter  responds  to  those  excitations,  and  how  those  responses  can 
be  changed.  Force  determination  divides  the  total  problem  into  its 
definable  and  solvable  parts,  and  gives  the  engineer  information  pertinent 
to  dynamic  evaluation. 

The  force  determination  approach  divorces  the  unknown  and  undefinable 
classical  dynamics  of  the  rotor  in  the  fixed  system  from  the  well-defined 
classical  dynamics  of  the  fuselage,  thereby  eliminating  the  situation  and 
the  problem  mentioned  above.  Periodic  coefficients  or  nonlinearities  in 
the  fixed  system  equations  of  the  rotor  are  Immaterial  in  force  determina- 
tion. The  testing  is  rigorously  correct  when  using  zero  hub  mass  because 
the  response  of  the  fuselage  to  rotor  forces  and  moments  is  determined  as 
the  sum  of  products  of  the  hub  loads  and  the  partial  derivatives  of  the 
fuselage  response  with  respect  to  the  hub  loads.  The  undefined  question 
of  "total  helicopter  natural  frequencies"  thus  does  not  enter  into  the 
matter.  It  is  replaced  by  more  practical  considerations  of  maximum  fuse- 
lage response  (roughly  "fuselage  natural  frequencies")  to  rotor  excita- 
tions and  maximum  actual  rotor  excitations  from  flight. 

Full  Ship  Ground  "Flight"  Testing 

Once  having  used  force  determination  to  find  the  magnitudes  and  phasings 
of  the  vibratory  forces  and  moments  acting  on  the  aircraft  in  all  the 
major  flight  conditions  of  the  mission  profile,  engineers  can  reproduce 
these  loads  on  an  actual  ship  in  a test  bay.  The  reliability  of  the 
engine  "motored"  (turning  slowly),  the  main  and  tail  rotor  transmission 
mounts  and  cases,  the  fuselage  structure,  the  stabilization  equipment, 
the  secondary  structure,  the  radios,  the  flight  instruments,  the  naviga- 
tion equipment,  the  external  stores  supports  and  other  parts  of  the 
aircraft  affected  by  the  main  and  tall  rotor  Induced  vibration  flight 
environment  can  then  be  evaluated  by  "ground  flying"  the  ship  24  hours  a 
day,  7 days  a week. 

Not  only  would  this  "ground  flying"  provide  a very  significant  calendar 
time  lead  over  actually  flying  aircraft  to  accumulate  fatigue  damage  data, 
at  the  rate  of  up  to  eight  thousand  seven  hundred  flight  hours  per  year, 
but,  because  the  ship  Is  on  the  ground,  conditions  pertinent  to  the 
damages  can  be  easily  observed.  The  opportunity  to  make  such  ground 
measurements  as  passive  vibration  surveys,  acoustic  emission  monitoring 
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and  electronic  measurements  with  the  ship  in  its  flight  vibration  environ- 
ment would  open  unprecedented  possibilities  for  determining  the  causes  of 
failures  and  for  testing  potential  fixes  long  before  the  failures  occur 
in  flight. 

The  advantages  of  "ground  flying"  reliability  testing  apply  to  ships  with 
a long  service  record  as  well  as  to  new  aircraft  under  development. 

Changes  in  mission  profile  and  minor  changes  in  aircraft  equipment  and 
structure  often  result  in  new  reliability  and  maintainability  problems, 
many  of  which  could  be  prevented  if  the  ship  could  be  externally  observed 
in  its  flight  vibration  environment.  The  kind  of  "ground  flying"  testing 
which  force  determination  makes  possible  is  by  no  means  all-inclusive. 

The  effects  of  main  and  tail  rotor  vibratory  forces  and  the  effects  of 
other  vibratory  loadings  on  most  of  the  helicopter  components  can  be 
simulated, but  ground  flying  does  not  include  flight  effects  on  the  rotors. 
Rotating  parts  of  the  transmission,  at  least  under  torque  load,  are  not 
accounted  for,  but  the  fatigue  effects  on  the  transmission  case,  which  is 
a major  load  path,  are  simulated.  The  effects  of  the  rotor-induced  vi- 
bration environment  on  the  actual  engine  can  be  ground-flight  tested 
presuming  that  the  engine  is  lubricated  and  motored,  perhaps  externally 
under  no  output  load,  to  prevent  false  brinelling.  Engines  on  helicopters 
have  been  shaken  with  the  engine  externally  motored.  In  "ground  flying", 
the  helicopter  would  be  under  the  steady  flight  1.0-g  download  as, 
initially,  the  limit  of  interaction  of  steady  and  vibratory  loadings. 
Eventually,  using  the  Elastic  Influence  Coefficients  to  accurately 
simulate  greater  than  1.0  g-steady  loads,  engineers  can  further  refine 
the  total  maneuver  load  simulation  in  "ground  flying". 


The  inclusion  of  all  main-  and  tail-rotor  harmonics  of  rotor  speed  in 
"ground  flying"  is  not  beyond  the  state  of  the  art,  but  insofar  as  simu- 
lating principal  effects  at  reasonable  cost,  it  should  perhaps  be  left 
to  eventual  refinement. 


How  to  "Ground  Fly"  a Helicopter 

"Ground  flying"  a helicopter,  even  given  the  vibratory  rotor  loads  from 
force  determination,  may  not  seem  to  be  an  entirely  straightforward  task, 
but  It  is  quite  within  the  state  of  the  art.  The  first  requirement  is 
that  we  reproduce  specific  shears,  moments  and  respective  phase  angles 
at  the  hub  using  shakers..  For  an  arbitrary  phase  reference,  say  "A" 
blade  passing  zero  azimuth,  we  have  twelve  variables  for  each  harmonic 
at  each  hub:  three  forces,  three  moments  and  a phase  angle  associated  with 
each.  At  the  hub,  this  condition  can  be  satisfied  by  six  shakers  properly 
arranged  and  properly  controlled.  The  technology  for  the  control  of  the 
phase  and  force  magnitude  of  an  electromagnetic  shaker  relative  to  any 
other  electromagnetic  shaker  In  a system  is  well  established  and  equip- 
ment for  doing  this  Is  commercially  available;  such  equipment  was  used 
in  this  contract.  The  essential  problem  in  "ground  flight"  testing  Is 
then  seen  to  be  that  of  using  the  mathematically-necessary  six  shakers 
In  a manner  that  causes  negligible  interaction  among  the  various  co- 
ordinates in  essentially  reproducing  the  shears,  moments  and  associated 
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phase  angles  found  in  force  determination  in  flight. 


Figure  2 shows  a method  of  doing  this.  Essential  to  this  particular 
method  are  relatively  long  cables  and  a lightweight  pyramid  which  trans- 
mits forces  to  the  hub.  The  long  cables  minimize  interaction,  and  the 
lightweight  three-dimensional  truss  construction  of  the  force-transmitting 
pyramid  minimizes  drainage  of  the  shaker  forces  to  inertial  space  while 
retaining  essential  strength.  The  cables  link  the  aircraft  to  the  shaker 
with  low  spring-rate  shock  cord  at  both  extremities  beyond  the  shaker 
armature/aircraft  connection.  The  armatures  of  the  shakers  are,  of 
course,  structurally  integral  with  the  cables.  The  cables  are  preloaded 
to  static  loads  in  excess  of  the  vibratory  forces  produced  by  the  shakers, 
thereby  allowing  the  shakers,  in  the  phraseology  of  freshman  engineering, 
to  "push  on  a rope".  The  low  spring  rate  of  the  shock  cord,  or  "bungee", 
extreme  to  the  aircraft  and  the  shakers,  makes  the  mechanical  impedance 
load  of  the  tare  on  the  shakers  negligible  compared  to  the  mechanical 
impedance  load  of  the  aircraft  on  the  shakers  for  all  practical  purposes. 
The  spring  rates  of  the  cables  do  not  diminish  force  transmission  and 
are  important  only  to  the  extent,  in  the  frequency  range  of  normal  heli- 
copter testing,  that  they  are  not  so  low  as  to  not  require  excesses  of 
the  travel  limitations  of  the  shaker  armatures.  This  matter  is  one  of 
designer  control,  and  solid  steel  straps  with  rod-end  bearing  connections 
might  serve  as  well  as  cables. 

The  advantage  of  long  cables,  or  straps,  in  minimizing  interaction  among 
the  coordinates  is  best  illustrated  by  example.  With  cables  20  feet  long 
on  either  side  of  the  pyramid  preloaded  statically  to  3000  pounds,  a hub 
vibration  of  1.0  q at  10  cps  produces  a transverse  interaction  load  of 
only  2.45  pounds.  For  a given  frequency  and  g-loading  at  the  hub,  the 
transverse  interaction  load  is  inversely  proportional  to  the  equal  cable 
length  on  either  side  of  the  ship.  One  would  preload  the  cables  to 
slightly  greater  than  the  maximum  vibratory  force  to  be  transmitted  to 
the  hub.  With  a 20-foot  half  span  in  a static  test  bay,  we  should  there- 
fore expect  at  1.0  g and  10  cps  an  interaction  of  something  less  than  one 
part  in  a thousand. 

In  a teetering  rotor  helicopter,  moments  other  than  torque  are  negligible, 
and  the  setup  for  "ground  flying"  is  somewhat  simpler,  as  shown  in 
Figure  3.  For  simplicity,  the  shakers,  shock  cord,  tensioners  and  other 
devices  have  been  eliminated  in  Figure  3. 

The  forces  and  phase  angles  for  each  shaker  necessary  to  produce  the 
desired  moment,  shear  and  associated  phase  angles  are  given  by  a simple 
4x4  matrix  inversion.  With  the  distances  x,  y,  and  z measured  from 
the  hub  center,  the  real  and  imaginary  forces  on  shakers  1 and  2 to 
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Figure  2.  Method  of  Ground  Flying  a Helicopter 
With  More  Than  Two  Blades 


Figure  3.  Method  of  Ground  Flying  a Two-Bladed  Teetering  Rotor  Helicopter 
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produce  the  pitching  moments  MR  and  M*  and  the  vertical  shears  FR  and  F1 
are  given  by:  y y 2 z 


By  way  of  illustration,  let  = x2  = 30  inches,  and  assume  the  pitching 

moment  to  be  10,050  inch-pounds  at  a phase  angle  of  5.7  degrees  and  the 
vertical  shear  to  be  2062  pounds  at  a phase  angle  of  14  degrees.  From 
Equation  1 it  is  found  that  shaker  1 must  produce  1197  pounds  at  a phase 
angle  of  15.6  degrees. 

An  equation  similar  to  Equation  1 can  be  written  by  inspection  for  each 
shear  and  moment  to  find  the  forces  and  phase  angles  of  each  pair  of 
shakers. 

Although  not  shown  in  Figures  2 and  3,  external  forces  at  points  other 
than  the  hub  can  be  readily  applied  in  the  "ground  flight".  Such  forces 
might  come  from  oscillatory  aerodynamic  loads,  vibratory  downwash 
impingement,  gunfire  or  other  external  loadings  that  can  be  found 
through  force  determination  and  may  be  included  in  the  "ground  flights" 
over  a mission  profile. 

Detecting  Unsuspected  Vibratory  Forces 

Force  determination  may  be  used  to  detect  and  measure  vibratory  excitation 
forces  from  flight  data  on  the  aircraft  that  may  not  have  been  anticipated 
as  disturbances.  In  an  actual  case  described  later  in  this  report,  the 
oscillating  angle  of  attack  of  the  horizontal  stabilizer  of  a UH-1H  was 
passing  a two-per-rev  vertical  excitation  from  the  swashplate  into  the 
tail  boom  during  a relatively  high  speed  forward  flight.  Needless  to 
say,  this  suggests  that  the  cockpit  vibration  of  such  helicopters  at  high 
speeds  might  be  reduced  simply  by  a change  in  the  control  interaction  of 
the  horizontal  stabilizer. 
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Force  determination  also  offers  the  opportunity  to  do  flight  investiga- 
tions of  the  vibration  effects  of  downwash  impingement  on  large  surfaces 
close  to  the  rotor,  the  effects  of  tail-rotor  fuselage  interference 
effects,  harmonic  forcing  from  aerodynamics  on  control  surfaces,  and  the 
effects  of  other  excitations  that  have  never  been  measured. 


PRINCIPLES  OF  FORCE  DETERMINATION 
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BASIC  APPROACH 


The  accelerations  of  a structure  at  any  frequency  are  given  by  the 
product  of  the  acceleration  mobility  matrix  and  the  external  forcds 
acting  on  the  structure. 

{y}  = [Y]{f } (2) 

The  accelerations  and  the  forces  are,  in  general,  complex;  i.e.,  they 
have  various  phase  angles,  and  therefore  the  acceleration  mobility  matrix 
is  also  complex.  As  used  here,  "forces"  should  be  taken  to  include 
moments  as  well  as  shears. 

If  there  are  N external  forces  and  moments  to  be  determined,  it  is 
necessary  to  use  at  least  N accelerometers  so  that  there  will  be,  in 
Equation  2,  as  many  simultaneous  equations  as  there  are  unknowns. 

These  accelerometers  will  not  be  at  the  points  of  application  of  the 
external  forces  but  will  be  distributed  throughout  the  fuselage.  The 
mobility  matrix  is  obtained  by  shake  testing  the  helicopter  fuselage 
without  the  rotor.  The  terms  in  the  mobility  matrix  are  the  real  and 
imaginary  ratios  (partial  derivatives)  of  the  response  of  each  accelerom- 
eter to  each  external  force  or  moment,  often  expressed  in  units  of 
g/pound.z  If  there  are  N forces  and  moments  to  be  found  in  flight  and 
N accelerometers  on  the  ship,  then  the  forces  and  moments  would  be  given 
by: 

(f}  = [Y]"1 {y}  (3) 


SOLUTION  TO  THE  SIGNAL/NOISE  PROBLEM 

Although  mathematically  correct  and  substantially  what  has  been  attempted 
with  hub  strain  gaging  in  the  past.  Equation  3 is  not  suitable  for 
practical  engineering  application  because  it  is  sensitive  to  errors  in 
the  measurement  of  accelerations  in  flight  and  of  mobilities  in  the 
ground  shake  test.  As  shown  in  Reference  6,  matrices  which  are  numeri- 
cally well -conditioned  may,  because  of  errors  in  the  measurement  of  their 
elements,  give  physical  information  which  is  meaningless  when  inverted. 


1 

.1 


■-i 


2 


Flannelly  and  Giansante. 


6 


W.G.  Flannelly,  A.  Berman  and  N.  Giansante,  RESEARCH  ON  STRUCTURAL 
DYNAMIC  TESTING  BY  IMPEDANCE  METHODS,  Kaman  Aerospace  Corporation; 
USAAMRDL  Technical  Report  72-63A,  B,  C,  & D,  U.  S.  Army  Air  Mobility 
Research  and  Development  Laboratory,  Fort  Eustls,  VA,  November  1972, 
AD  756389,  AD  756390,  AD  756391,  & AD  756392. 
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Further,  one  would  have  no  way  of  determining  from  the  use  of  Equation  3 
whether  the  information  obtained  is  physically  meaningless. 6 

To  eliminate  this  serious  problem,  we  use  the  property  of  "redundancy", 
which  has  been  employed  so  effectively  in  information  theory  and  has  been 
described  as  that  property  which  "...makes  language  intelligible  in  the 
midst  of  noise". o To  achieve  a high  degree  of  redundancy  we  should  like 
to  use  many  more  accelerometers  in  the  fuselage  than  there  are  forces  to 
be  determined.  This  would  make  the  mobility  matrix  of  Equation  2 
rectangular,  with  many  more  rows  than  columns,  giving  us  many  more  equa- 
tions than  unknowns.  The  problem  now  is  how  to  invert  a rectangular 
matri x. 

In  1968,  A.  Berman  of  Kaman  Aerospace  Corporation  discovered  a matrix 
form  with  peculiar  properties.  It  was  later  found  that  the  same  proper- 
ties could  be  achieved  with  a process  of  matrix  algebra  known  as  the 
"pseudoinverse",  "Moore's  generalized  inverse",  or  the  "generalized 
reciprocal".  Among  the  simplest  properties  of  the  pseudoinverse, 
discovered  by  E.  H.  Moore  in  1920  but  lost  and  rediscovered  by  R.  Penrose 
in  1955,  is  that  the  inversion  of  a rectangular  matrix  of  maximal  rank 
amounts  to  a least-squares  solution  to  a problem  with  more  equations  than 
unknowns. l 0 The  pseudoinverse  was  employed  in  dynamics  in  References 
2,  6,  and  11,  and  additionally  in  eigenvalue  form  for  the  first  time  in 
Reference  5,  where  a derivation  is  shown. 


d Flannel ly  and  Giansante. 
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Cambridge  Philosophical  Society.  51,  1955,  pp  406~4T3. 

^A.L.  Klosterman  and  J.R.  Lemon,  DYNAMIC  DESIGN  ANALYSIS  VIA  THE 
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If  we  now  use  F accelerometers  in  the  fuselage,  where  the  number  of 
accelerometers  used  is  much  larger  than  the  H forces  and  moments  at 

“ 1 

the  hub  (F  » H),  then  the  square  [Y]pH  inverse  matrix  of  Equation  3 

becomes  a rectangular  H x F pseudoinverse  matrix  [Y]pH*  The  superscript 

"+"  replaces  the  familiar  "-1"  to  denote  the  pseudoinverse.  We  thus 
obviate  the  problem  of  inversion  leading  to  meaningless  information, 
and  the  signal-to-noise  problem  is  solved.  Then  we  have,  from 
Equation  3: 


{F}h  = 

Hxl  HxF  Fxl 
(F»H) 


WRONG  SIGNALS  AND  NOISE 


The  usefulness  of  the  information  from  dynamic  testing,  or  analysis,  in 
the  process  of  deciding  how  to  improve  a helicopter  is  obviously  heavily 
dependent  on  the  physical  accuracy  of  the  information.  Hardware 
decisions  based  on  false  information  from  tests  and  analyses  which 
cannot  be  or  are  not  verified,  and  the  large  number  of  pointless  dynamic 
tests  conducted  annually  drain  the  limited  budget  of  funds  which  could 
be  more  usefully  applied  elsewhere. 

Although  the  signal-to-noise  problem  is  fundamental  both  to  dynamic 
testing  and  communication  theory,  dynamic  testing  differs  from  cormuni ca- 
tion in  that  it  involves  the  additional  possibility  of  receiving  "wrong 
signals".  Wrong  signals  can  come  from  many  common  causes:  a mounting 
distortion  of  an  accelerometer  case,  an  accelerometer  mounted  upside 
down  with  respect  to  the  physical  meaning  given  to  its  signal,  incorrect 
calibrations,  and  location  errors.  Electronic  or  statistical  checks 
against  noise  cannot  insure  that  a test  signal  is  physically  meaningful. 
The  only  way  a wrong  signal  can  be  found  is  through  checking  the  physics 
of  the  structure  itself.  Two  ways  to  do  this  in  helicopter  dynamic 
testing  are  Maxwell's  Reciprocity  Theorem  and  the  Rigid  Body  Acceleration 
Coefficient. 
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METHODS  OF  OBTAINING  THE  MOBILITY  MATRIX 


DIRECT  SHAKING 

The  most  direct  method  of  obtaining  the  mobility  matrix  is  to  shake  the 
helicopter  at  the  locations  of  the  excitation  forces  and  measure  the 
resulting  accelerations  per  unit  shaking  force  in  the  fuselage.  In  the 
determination  of  rotor  forces  and  moments,  this  usually  entails  shaking 
at  the  hub  with  vertical,  lateral,  and  longitudinal  forces  and  moments 
one  at  a time. 

Among  the  methods  discussed  in  this  report,  direct  shaking  has  the 
advantages  of:  providing  the  mobility  matrix  directly  without  any 
intervening  techniques;  exciting  the  fuselage  at  the  sources  of  flight 
excitation,  and  thus  allowing  shake  testing  consistent  with  known  dynamic 
testing  practice.  It  has  the  disadvantages  of:  the  customary  expensive 
set-up  costs;  the  problems  and  expense  of  minimizing  tare  dynamics  in 
moment  shaking;  the  necessity  to  accurately  cross-control  two  shakers 
simul taneously  to  simulate  moments;  and  the  lack  of  opportunity  to  use 
smaller-than-usual  shakers. 

Although  applicable  to  all  types  of  helicopters,  direct  shaking  is  most 
easily  applied  to  teetering  rotor  ships  where  the  only  significant 
vibration  hub  excitations  are  vertical,  fore  and  aft,  lateral  shears  and 
torque. 

RECIPROCITY  SHAKING 


This  is  not  only  another  method  of  obtaining  the  mobility  matrix,  but 
can  and  should  be  used  to  a limited  extent  with  other  methods  as  a 

physical  check  on  the  mobilities  measured.  The  basis  of  this  technique 

is  the  well-known  reciprocity  theorem  of  Maxwell.  It  states  that  in  any 
linear  system  the  acceleration  at  j due  to  a force  at  k (and  only  at  k) 
is  the  same  as  the  acceleration  at  k due  to  the  same  force  at  j.  In 
mobility  terms,  If  the  reciprocity  theorem  is  not  confirmed 

by  shaking  at  some  accelerometer  station,  the  engineer  should  carefully 
reexamine  his  work  for  either  an  error  in  the  testing  or  structural 

nonlinearity,  such  as  ankylosis,  which  can  and  should  be  accounted  for. 


R.A.  Frazer,  W.J.  Duncan,  and  A.R.  Collar,  ELEMENTARY  MATRICES, 
Cambridge,  England,  Syndics  of  the  Cambridge  University  Press, 
1938,  reprinted  1965. 


To  obtain  the  mobility  matrix  by  reciprocity  shaking,  the  engineer  would 
shake  at  each  of  the  accelerometer  stations  and  measure  the  responses 
(using  rotational  as  well  as  translational  accelerations)  at  the  excita- 
tion's coordinates.  The  advantages  of  this  method  are  that  the  suspension 
system  is  greatly  simplified  because  there  is  no  hub  shaking  and  smaller 
shakers  can  be  used  if  the  fuselage  accelerometers  are  placed  at  mobility 
antinodes  (places  of  high  fuselage  response).  The  main  disadvantage  of 
this  method  is  that  it  requires  shaking  at  the  many  accelerometer  loca- 
tions. This  requires  not  only  much  more  testing  time  than  direct  shaking 
but,  by  presently  proven  methods,  also  requires  large  set-up  costs  in 
providing  the  many  shaker  hard-point  attachments. 

At  this  time  it  is  difficult  to  envision  reciprocity  testing  having  an 
overall  advantage  as  a primary  method  for  obtaining  the  mobility  matrix 
over  the  other  two  methods  presented  here.  However,  impact  testing  with 
Fast  Fourier  Transform  equipment  has  made  very  significant  strides  in 
recent  years.  Should  this  technique  be  proven  practical  for  helicopters, 
then  reciprocity  testing  could  become  a very  economical  method  for  ob- 
taining the  mobility  matrix.  In  the  meantime,  the  limited  use  of 
reciprocity  testing  as  a physical  check  on  other  methods  is  very  valuable. 

MODAL  ACCELERATION  TESTING 


Modal  acceleration  testing  (MAT)  can  define  the  total  N x N mobility 
matrix,  for  N accelerometers,  by  shaking  at  only  one  point  on  the 
fuselage.  This  allows  force  determination  to  be  done  with  only  one 
shake  test.  The  degree  of  flexibility  inherent  in  the  one-shake- test 
calibration  allows  the  engineer  to  evaluate  proposed  modifications,  such 
as  the  addition  of  missiles  or  other  external  stores  and  changes  in 
rotors  and  rotor  speeds,  without  repeating  the  shake  test.  MAT  also 
yields  the  true  undamped  natural  frequencies  (as  opposed  to  mobility 
peaks)  and  the  actual  normal  modes  (as  opposed  to  mere  mobility  ratios). 
These  two  results  are  necessary  for  the  comparison  of  finite  element 
dynamics  predictions  with  test  results.  MAT  also  provides  the  basic 
information  necessary  for  system  identification. 2. 6 


In  actual  practice,  MAT  is  used  at  several  excitation  points  for  a phys- 
ical check  on  the  test  data  and  to  take  advantage  of  modal  density 
reduction  caused  by  near  symmetry  of  the  fuselage.  Because  the  excitation 
can  be  at  a mobility  antinode,  a minimum  capacity  shaker  can  be  used. 
Because  there  need  be  only  one  to  three  excitation  points,  which  may  be 
conveniently  chosen,  an  elaborate  suspension  structure  or  a hangar  static 
test  bay  is  not  necessary.  Consequently,  very  little  attachment  hardware 
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need  be  designed  or  built.  The  entire  set-up  may  consist  of  a fork  lift 
and  an  A-frame  with  bungee  connected  to  the  hub.  The  setup  costs  of 
MAT  are  quite  small,  and  the  operational  costs  are  minimal. 

Because  of  the  importance  of  modal  acceleration  testing,  it  will  be 
described  in  detail  in  separate  sections  of  this  report. 


r 1 
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MODAL  ACCELERATION  TESTING  THEORY 


MODAL  ACCELERATIONS 


Derivation  of  Basic  Equation 

The  theory  and  application  of  analytical  testing  and  system  identification 
technology  include  many  mathematical  formulations  which  describe  the 
dynamic  behavior  of  linear  structures  having  proportional  structural 
damping.*3  At  any  frequency  w,  the  resulting  equations  of  motion  can 
be  written  as: 


[K  - u>  M + igK]  {y}  = {f} 


(5) 


The  terms  in  the  matrix  of  Equation  5 are  defined  as  the  "displacement 

C •• 

impedance"  Z.  Similarly,  the  "acceleration  impedance"  Z is  defined  as 
the  matrix  coefficients  of  the  accelerations  in  the  equations  giving 
the  forces 


[Z]  (y } = {f} 


(6) 


Note  that  Z is  not  the  second  derivative  of  Z with  respect  to  time,  which 
is  self-evident  in  the  context  it  is  used. 


From  Equation  6,  the  accelerations  are  given  by 

(y>  = [z]_1{f}  = [Y ] { f 1 


(7) 


P'-l 


where  the  inverse  of  Z is  defined  as  the  "acceleration  mobility"  Y.  Un- 
less otherwise  specified,  the  use  of  the  term  "mobility"  is  taken  to  mean 
"acceleration  mobility".  Helicopters  are  generally  shake-tested  while 
suspended  as  free  bodies,  and  "acceleration  mobility"  is  finite  at  zero 
frequency.  For  this  reason  and  others  which  become  obvious  from  the 
following  development,  "acceleration  mobility"  is  chosen  as  the  most 
convenient  "measurable"  quantity. 


3 Cransdale,  Gaukroger,  and  Skingle. 

13W.W.  Soroka,  NOTE  ON  THE  RELATIONS  BETWEEN  VISCOUS  AND  STRUCTURAL 
DAMPING  COEFFICIENTS,  Journal  of  the  Aeronautical  Sciences,  Volume 
16,  No.  7,  July  1949,  p"109; 
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As  seen  from  the  matrix  expressions  of  Equations  6 and  7,  both  impedance 
and  mobility  are  properly  partial  derivatives.  Acceleration  impedance  is 
the  partial  derivative  of  force  with  respect  to  acceleration,  and 
acceleration  mobility  is  the  partial  derivative  of  acceleration  with 
respect  to  force. 

The  acceleration  impedance  is  written: 


Z = M - (1  + ig)K/</ 

T T 1 

Premultiplying  Z by  $ 4>  and  postmul tiplying  by  $ $ gives 
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and  Ajxi  1S.  defined  as  the  Jx-th  modal  acceleration  of  the  i-th  mode. 


Letting 


v.  . - 

Jxl  ' «>{[»]{*>, 


J o - 

!1^[0  - 0)^/£j|)Z  + g?] 


id)  7? 


nf  [0  - w‘  vnfr  + gj] 


Equations  12  and  13  may  be  written  more  concisely  as 


Z Ajxlffc 

1=1 


••  T N 

YjXO.  = E AjxiF{c 


The  modal  acceleration,  Ajxj>  is  a term  which  can  be  determined  from  the 

plot  of  mobility  versus  frequency  for  the  acceleration  at  j due  to  a force 
at  x. 

It  is  more  convenient  to  use  the  orthonormal  mode  <p,  defined  as 

GpjJCMKi^  * [I]  (17) 
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than  the  conventionally  normalized  normal  mode  $>.  For  the  i-th  mode, 
then, 

Ajxi  = ^xt^jl  (,8> 


For  the  accelerometer  at  the  excitation  point  x,  we  determine  Axx-  for 
each  mode  i and  find 


For  any  accelerometer  j. 


\p  . = Ja 

YX1  > XXI 


09) 


rJi 


J A 7 

If  XXI 


(20) 


Suppose  that  the  shaker  is  at  station  x in  the  shake  tests  and  that  we 
wish  to  know  the  response  at  accelerometer  j if  the  ship  were  shaken 
at  accelerometer  k.  From  the  plot,  we  find  the  modal  accelerations 

Axxl-  AxxII-  AxxIII’  — ! From  the  Yjx  p,ot'  “ f1nd  AjxP  AjxII* 

AjxIII’  •••’  and  fr0ln  the  Vkx  plot’  “e  f,nd  Akxl’  Akxir  Akxl I I 

Using  Equations  19  and  20,  we  determine  the  orthonormal  mode  elements 
^xr  ^xir  ^xiir  •••*  *ji»  ^jir  ^jiir  •••  and  ^kir  ^kiir*** 
Then  we  construct 


YR  + iY1  = 
jkcj  lYjk 


> -T 

E Ajki<FL + tFl> 

i=l 


"n  "t 

E Vki(FL + iFl> 

i=l 


(21) 


Using  Equation  20,  the  engineer  can  construct  the  mobility  relative  to 

any  two  instrumented  points.  If  we  were  determining  the  acceleration  at 

•• 

j due  to  a moment  at  k in  flight,  then  the  plot  Ykx  of  the  shake  test 

would  be  that  of  the  rotational  acceleration  at  k due  to  a shaker  force 
at  x. 


ft 


As  shown  in  References  2 and  6,  and  as  can  be  seen  from  Equations  12 
and  13,  the  higher  the  natural  frequency  of  any  mode  above  the  fre- 
quency range  of  interest,  the  less  effect  it  has  on  the  mobility.  There- 
fore, the  number  of  modes  which  must  be  covered  in  the  testing  is  rela- 
tively small.  Usually,  covering  only  two  modes  above  the  frequency  of 
the  force  to  be  determined  (N-per-rev  in  most  rotor  cases)  is  sufficient. 

In  an  unrestrained  system,  such  as  a helicopter  in  flight,  there  are  six 
"rigid-body"  degrees  of  freedom.  Some  systems  are  not  totally  unre- 
strained, such  as  a beam  with  one  pin  end,  and  may  have  fewer  than  six, 
say  z,  "rigid-body"  modes.  Consistent  with  Reference  2,  Equation  15 
may  be  written 


YR  = E . 
jxoj  jx 


N 

£ 

i=Z+l 


A — 

j*1  a? 


2/02 

u)  /n. 


d - wvajr  + gf 


where  E.  is  the  rigid-body  acceleration  coefficient  (RAC).  At  zero 

j X 

forcing  frequency,  the  real  acceleration  mobility  is  equal  to  the  RAC. 
The  RAC  has  no  effect  on  the  imaginary  mobility.  The  RAC  can  be 
accurately  determined  from  weights  data  on  the  aircraft,  and  one  of  the 
checks  on  the  physical  validity  of  the  mobility  test  data  is  to  see 
whether  the  mobility  asymptotically  approaches  the  RAC  at  frequencies 
between  the  very  low  suspension  modes  and  the  first  elastic  mode  of  the 
fuselage. 

Anti  resonance  Theory  Representation 

The  following  development  is  a completely  mathematical  formulation  of 
Yjx  based  entirely  upon  antiresonant  eigenvalues  of  matrices.^ 
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Mathematically,  the  jx-th  acceleration  mobility  at  frequency  id  is  defined 
as 

ay,-  H)j+X|z|ix 

Y = — J...  = A* 

OX  3FX  |?|  (23) 

since  Z"^  = Y where  |Z|  is  the  determinant  of  Z and  (-1)^+X|Z|  . is  the 

••  J X 

jx-th  co-factor  of  Z. 

The  acceleration  impedance  from  Equation  7 can  be  written  as 


Z = m[  [I]  + M-1  (1  + ig)K/u)2  ] 


The  determinant  of  Z is  then  seen  to  be 


|Z|=  | M | | [I]  + M-1 (1  + ig)KA>2)| 


dhich,  when  expanded,  becomes 


|Z|  . a,2N  - p,*2*"'1' 


. - (-i )N_1  n + m 


vhere  N is  the  order  of  the  matrix  and 


|K|  N 2 

— = n 

|M|  " 1 

3^  are  the  eigenvalues  of  M"^K  and  Equation  26  is  further  simplified  as 


|z|  » -14  TT  of  [0  ♦ 19)  - “2/«il 

id  i*l  L J 


Similarly,  the  jx-th  co-factor  of  Z becomes 


H)j+Xizijx  = ^f  naj2xi  [ci  * 19)  - «2/«jx1]  (28) 
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where  a^x  are  the  eigenvalues  of  l^ljxl^ljx  known  as  anti  resonances. 

Reference  14  interprets  the  meaning  of  a both  mathematically  and 
physically.  JX1 

In  an  unrestrained  system,  there  are  Z natural  frequencies  (ft.. ) which  are 
zero,  and  Z-l  antiresonant  frequencies  (ajx)  which  are  zero.  Combining 
Equations  27  and  28  gives 


From  Equation  29  the  mobility  can  be  reproduced  if  E.  , a..,  ft.,  and 
g are  known.  JX  J 

REMOVAL  OF  DAMPING 

By  definition,  the  natural  frequencies  and  modal  accelerations  are  un- 
damped parameters  of  the  system.  In  order  to  calculate  these  undamped 
modal  parameters,  the  undamped  mobility  must  be  known.  Since  all  measure- 
ments include  damped  real  mobility,  it  is  necessary  to  derive  the 
appropriate  relations  for  calculating  undamped  resonances.  As  a con- 
sequence of  the  development,  a new  expression  for  determining  modal 
damping  is  included.  This  new  result  is  valid  for  closely  coupled  modes 
2 

and  assumes  g^  « g^ . 

Consider  Equation  29  and  rewrite  as: 


I 


1 
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where  the  j,  k and  x subscripts  are  left  off  ai  for  clarity. 

Reindexing  the  anti  resonances  to  pair  with  the  resonances.  Equation  31 
becomes 

V Ejx  TT  (Ri  + iI1-)  , 

i=Z+l  1 

Ignoring  terms  involving  the  products  of  the  modal  imaginaries  in  the 
expansion  of  Equation  31,  the  real  mobility  can  be  approximated  as 


r<i  - v*1  ■ v + flfi 

yr  - e TT  di  ni  -I 

Jx  Jx  W+1  [0  - Jn*)*  * 4]  c 

where  the  j,  k subscripts  are  left  off  a^  for  clarity. 

Expressing  the  numerator  and  denominator  of  the  factors  of  Equation  33 
in  terms  of  the  roots  given, 

(*7  - “f)^  - bJ) 

~r  N fif  1 n.  1 

’Sx  ! Ejx  n jt1 yp 

% - (l  + ig)  \ - (1  - ig) 

ni  <3 


where 


If,  for  any  mode  i,  the  a.,  and  roots  are  complex,  then  there  will  not 

be  any  zero  mobility  crossover  for  that  mode.  To  have  a zero  crossover 
in  the  real  mobility  for  the  i-th  mode,  it  is  necessary  and  sufficient 
that,  for  a real  antiresonance,  the  quantity  under  the  radical  of  and 

8^  be  positive.  This  condition  gives  the  range  of  a^/ft..  in  terms  of  g. 

for  which  a zero  crossover  in  the  real  acceleration  mobility  exists. 

In  other  words, 


i + 9}  + 9i  ±n.  < H - gi  (37) 

When  there  is  a zero  crossover  for  a real  anti  resonance,  the  crossover  at 
ujg=  B^  is  associated  with,  but  not  equal  to,  the  resonance,  and  the 

crossover  at  u>a  = a.fij  is  associated  with,  but  not  equal  to,  the  anti- 
resonance. An  accurate  approximation  of  the  true  undamped  antiresonance 
is  obtained  by  adding  Equations  35  and  36. 


2-2  2 2 

ai  = \ + \ " n1  (38) 

Subtracting  Equations  35  and  36  gives  an  approximation  for  modal  damping: 


n2  - 
9i  ’ 

ni  % 


2 2 
CD 

6i  “i 


Equations  38  and  39  assume  that  ft..  is  the  undamped  resonance. 

The  determination  of  the  true  undamped  resonance  can  be  performed  as 
follows.  If  a constant  K is  subtracted  from  the  real  acceleration 
mobility,  Equation  33  becomes,  after  simplification. 


'jV'*  (Elk  - K)  ft  -L 

J 1-2+1  n? 


(CD  - CD  )((D  - (Do  ) 
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(1  - T>  *9? 
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where  cd  and  cdq  are  frequencies  at  which  the  mobility  is  equal  to  K. 
a1  Pi 

Let  cdu  be  the  frequency  at  which  the  mobility  equals  K for  no  damping; 
Ki 

then  the  roots  may  be  expressed  as 
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Adding  Equations  41  and  42  gives  the  frequency  at  which  the  mobility 
is  equal  to  K when  the  damping  is  zero. 


U)^  =(/+  U>1  - «? 

ki  “1  6i  ’ 


Similarly,  subtracting  the  two  equations  gives  the  damping  coefficient 

uj2  ui2 

.2  - 8i  . 

1 «?  <<,  * 4t  - "1  > (44) 

Equations  43  and  44  are  independent  of  K,  which  means  that  they  are  valid 
for  and  independent  of  the  mobility  at  which  the  damped  real  mobility 
exists  at  adjacent  frequencies.  The  calculation  of  the  undamped  mobility 
is  based  on  a frequency  shift  of  the  damped  mobility  and  does  not  yield 
undamped  mobilities  at  levels  where  the  damped  mobility  does  not  exist. 
Equations  41  and  42  can  now  be  used  to  determine  the  undamped  resonances. 
At  a maximum  or  a minimum  in  the  damped  real  mobility.  Equations  41  and 
42  are  equal . 
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Figure  4.  Real  Mobility 


This  occurs  when  the  term  under  the  radical  is  zero. 


(l  + wj;  /a\)  = 4(1  + g?)  c ^ /n? 


(46 


There  are  two  such  frequencies  for  each  mode  (p,  and  p2  ) in  the  damped 

i 

real  mobility,  as  shown  in  Figure  4.  From  Equation  43  this  gives 
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Equations  47  and  48  are  substituted  into  Equation  46  to  obtain 


Dividing  Equation  49  by  Equation  50  gives,  after  simplifying, 


2pl 

n2  = -A 

“l  n2 

V 


Equation  51  relates  the  true  undamped  resonance  in  terms  of  the  peak 
frequencies  p,  and  p~  independent  of  damping  and  mobility, 
i i 

Substituting  Equation  51  into  Equation  44  similarly  relates  the 
damping  solely  in  terms  of  these  modal  peak  frequencies. 


. (Pjjj/y2  - 1 , /p?<  p„\ 

1 2 p21/p„  - 2 ^p„.  - 
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Equations  51  and  52  are  both  important  in  the  applications  of  the  E-MASS 
for  reducing  the  data  in  terms  of  undamped  modal  accelerations. 
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The  following  derivation,  although  shown  in  Reference  2,  is  presented 
here  in  the  interest  of  completeness. 


In  Figure  5, x,  y,  and  z are  any  orthogonal  axes  through  the  center  of 
gravity  of  a body  of  mass  M,  with  points  x^  y and  z^,  shown  on  the 

respective  axes,  having  rotational  accelerations  0 , 0 , and  0 . 

x y z 


2 


Figure  5.  Axis  System 

The  moments  M , M , and  M about  the  x,  y,  and  z axes  may  be  written 
x y 4 

using  Lagrange's  Equation  as 
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where  Ix,  Iy,  and  Iz  are  the  moments  of  inertia  about  the  x,  y,  and  z 
axes,  and  U^,  Uxz,  and  U are  the  products  of  inertia  with  respect 
to  the  x and  y axes,  the  x and  z axes,  and  the  y and  z axes. 


Equation  53  may  be  written  in  general  notation  as 


{Mq} 


3M 
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The  [3M/30]  matrix  may  be  inverted  to  yield 
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(55) 


where  the  determinant  of  the  matrix  is  given  by  A = I I I -2U  U U 
p?2  x y z xy  zy  xz 

IT  I - IT  I - U I 
xy  z zy  x xz  y 


Solving  for  the  rotational  accelerations  and  using  the  definition  pre- 
sented in  Equation  55  yields 


(0)  = 


|{Mq} 


(56) 


The  matrix  [30/aM]  represents  the  rigid-body  acceleration  coefficients  of 
the  rotational  acceleration  about  any  of  the  specified  axes  with  respect 
to  a moment  about  any  of  the  given  axes. 

The  linear  accelerations  of  the  points  xi , y. , and  z..  due  solely  to 
rotations  or  the  rigid  body  may  be  expressed  as 


Equation  57  may  be  written  more  generally  using  Equation  56, 


{q>  = [Q^te)  = [Q^Cae/aMltM  } 


(58) 


The  moments  about  the  x,  y,  and  z axes  caused  by  orthogonal  forces  at  any 
point  are  given  by 
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(59) 


Equation  59  may  be  written  as 


(Mq)  = [Qp]{F} 


(60) 


Newton's  Law  relates  the  accelerations  of  the  center  of  gravity  of  the 
body  to  the  applied  forces  and  the  mass  of  the  body: 


<V  ■ »., 
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(61) 


The  rigid-body  acceleration  coefficients  describing  rotational  accelera 
tion  due  to  moment  are  given  by  Equation  55  and  are  written  in  matrix 
form  as 


^E9M^  = 


36 

3M 


The  rigid-body  acceleration  coefficients  relating  translational  accelera- 
tion and  moment  are  obtained  from  Equation  58: 


[Ejjm]  = [3q/3M]  = [Qi][30/3M]  = [Q.][EqM] 
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The  rigid-body  acceleration  coefficients  presenting  rotational  accelera- 
tion due  to  force  are  defined  from  Equations  56  and  60: 


CEqF]  = [36/aF]  - U^J  [Qp]  - [EgM][QF]  if* 

The  accelerations  of  a point  on  the  rigid  body  result  from  the  trans- 
lational acceleration  of  the  center  of  gravity  and  the  linear  accelera- 
tion due  to  pure  rotations  of  the  rigid  body.  The  translational 
acceleration  is  obtained  from  Equations  58,  60,  and  61  and  is 
written  as 


fq)  - tF)  Jf  + [Q^tE^HQpKF) 


Equation  65  yields  the  rigid-body  acceleration  coefficients  relating 
translational  acceleration  and  force: 


Pqf]  5 t3q73F]  > [Q,][E5m][Qf] 


* METHOD  FOR  THE  EXPERIMENTAL  DETERMINATION  OF 
THE  RIGID  BODY  ACCELERATION  COEFFICIENTS 

In  some  cases,  where  a commercial  helicopter  has  been  adopted  for  military 
use  or  where  a military  helicopter  has  undergone  extensive  modification, 
knowledge  of  the  principal  inertial  axes  and  moments  of  inertia  may  be 
doubtful.  If  we  determined  the  inertial  properties  (the  gross  weights 
and  center  of  gravity  location  are  known  from  weight  and  balance  measure- 
ment) by  swinqing  the  ship,  we  could  calculate  the  rigid  body  acceleration 
coefficients  (RACs)  by  Equation  66,  but  accurately  determining  the  in- 
ertial properties  by  swinging  a fuselage  is  extremely  difficult.  We  would 
determine  each  RAC  from  the  low  frequency  response  of  each  mobility  run, 
or  by  the  matrix  inversion  method  of  determining  modal  accelerations^ , 


R R 

or  by  subtracting  fr0m  Yw-j  (Equation  22).  However,  these  would 

not,  In  themselves,  give  us  Information  about  the  RACs  for  stations  at 
which  the  helicopter  was  not  shaken  in  the  ground  test.  But,  if  we  use 
the  method  described  below,  we  can  experimentally  determine  all  the  RACs 
by  shaking  at  only  three  points.  Since  the  RACs  derived  from  this  test 
would  be  independent  of  those  measured  in  the  shake  testing,  they  would 
serve  as  a physical  property  cross-check  on  that  data.  The  method  makes 
no  assumptions  about  symmetry  or  principal  axes. 


Klosterman  and  Lemon. 
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The  engineer  chooses  three  points  on  the  aircraft;  x^,  y1 , z1 ; x2,  y2, 

z2;  and  Xy  y3,  z3;  measured  from  the  center  of  gravity,  with  at  least 

y3>  z2,  and  y^  not  equal  to  zero.  By  any  of  the  methods  mentioned  above, 

the  engineer  finds:  E , E , E , E , E_  , and  E,  _ . Let  5 

11  21  Z3X1  X2X2  Z3X2  Z3Z3 

be  the  matrix  of  the  partial  derivatives  of  rotational  accelerations  with 
respect  to  applied  moments,  a 3 x 3 matrix.  Then, 


[E]  E [36/3M] 


It  is  seen  that 


where  M is  the  total  mass  of  the  ship 


We  therefore  see  that 
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The  three-dimensional, translational  RAC's  between  points  j and  k are  seen 
to  be  given  by 


w 


"°  zj  -yj“ 

-zj  0 xj 

[Q]-'  ([£•]-  [jr])[Q]-T 

■°  -zk  yk 
2k  0 -xk 

1 

o 

X 

1 

•r-> 

>> 

1 

_-yk  xk  0 _ 

xjxh 


yjxk 


zjxk 


'xjyk 

yjyk 


zjyk 


'xjzk 


yjzk 


ZjZk 


(70a) 


The  rotational/moment  RAC's  are  given  by  Equation  69.  Equation  70a 
shows  the  common  translation'force  RACs.  The  translation/moment  and 
rotation/ force  RACs  may  be  obtained  from  Equations  63,  64,  and  70 
as  follows: 

V - BMM 


[EjjF]  = [5][Qf] 


(70b) 

(70c) 


A detailed  description  of  the  pseudoinverse  iteration  method  for  obtaining 
modal  accelerations  and  an  experimental  verification  of  the  method  are 
given  in  Reference  2.  Other  methods  using  matrices  may  be  found  in 
Reference  6 and  Reference  11.  More  recently  developed  methods,  which 
are  considerably  simpler  to  use,  are  discussed  below. 

ANTI  RESONANCE  THEORY  METHODS 


Driving  Point  Check  Method 


The  following  develops  the  mathematical  expression  for  modal  acceleration 
which  was  presented  without  proof  in  Reference  2. 

The  proof  is  based  upon  antiresonance  theory,  and  it  is  shown  that  modal 
acceleration  can  be  determined  if  the  resonances,  the  anti  resonances , and 
the  rigid-body  acceleration  coefficients  are  known. 

In  Equation  22,  for  a force  at  k,  instead  of  at  x,  let  g = 0,  which  is 
the  case  for  an  undamped  system.  This  yields 
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where  the  superscript  "R",  indicating  the  real  part,  is  dropped  because 
there  is  no  imaginary  mobility  in  an  undamped  system.  Similarly,  for  an 
undamped  system  (g  = 0),  Equation  29  becomes 
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Flannel ly  and  Giansante. 

Flannel ly,  Berman,  and  Giansante. 


Klosterman  and  Lemon. 


Equation  77  may  be  written  in  the  form, 
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Equating  the  two  expressions  for  Yjk  given  by  72  and  73, 
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When  a)  = where  m is  one  of  the  elastic  modes  and  (Z+l)  £ M _<  N, 
Equation  74  becomes 
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Solving  for  Ajk^  gives  the  expression  reported  without  proof  in  Reference 
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The  anti  resonant  frequencies  of  a driving  point  are  all  real,  and  this 
equation  may  be  used  directly  as  an  independent,  hand-calculator  check 
on  the  driving  point  modal  accelerations  whenever  the  damping  is  low 
enough  to  detect  the  approximate  anti  resonant  frequencies. 
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Vector  Angle  Method  for  Transfer  and  Driving  Points 

Now  we  will  show  that  the  undamped  mobility  may  be  found  using  the  rigid 
body  acceleration  coefficient,  the  natural  frequencies,  and  the  anti- 
resonant frequencies.  However,  the  anti  resonant  frequencies  are  not 
necessarily  all  real  in  a transfer  mobility  plot;  there  may  be  anti- 
resonant roots  in  the  imaginary  frequency  spectrum.  Equation  76  cannot 
be  solved  without  all  the  antiresonant  frequencies,  and  imaginary  anti- 
resonances cannot  be  found  in  testing.  It  will  be  shown  that  undamped 
mobilities  may  be  obtained  using  only  the  real  anti  resonances. 

In  the  interest  of  notational  simplicity,  the  elastic  modes  are  numbered 
starting  with  "1",  instead  of  Z+l,  in  Equation  29. 

The  damped  mobility,  with  structural  damping,  in  a free  body  is  given 
exactly  by  Equation  29  for  an  infinite  number  of  modes: 
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Each  term  in  the  numerator  may  be  considered  to  be  a vector  of  length 
|(1  - o)2/ajki ) + g?|  and  angle  viu)  = arctanjg./(l  - tA>2/a2ki  )J  . 
Similarly,  the  denominator  is  the  product  of  vectors  of  length 

.2  ,r2\ 
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|(1  - oj2/Q2 ) + gZ|  and  angle  6.^  = arctan|g../(l  - wZ/ftZ)j 
Equation  29  may  be  expressed  as 
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Yjk  = Ejk 
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Note  that  the  absolute  amplitude  of  the  mobility  in  terms  of  the  sign  of 
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where  SGN  x = + for  x > 0 
= - for  x < 0 


The  undamped  mobility  is 
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The  cosine  terms  concerned  with  anti  resonances  are 
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If  the  antiresonance  is  imaginary,  then  o)  /a^  is  always  positive  and  the 
value  of  cos  v,  = +1/,/,  . _2,  which,  considering  that  g is  always  much 


less  than  one,  is  always  approximately  unity.  The  minimum  value  of  cos  v 
is  approached  as  the  anti  resonant  frequency  becomes  very  large  compared  to 
the  forcing  frequency  in  the  case  of  an  imaginary  or  real  anti  resonance, 
but  this  is  also  seen  to  be  the  case  when  cos  6 is  in  the  denominator. 
Because  there  is  necessarily  the  same  number  of  real  or  imaginary  anti- 
resonances  in  a free  system  as  natural  frequencies  (one  less  in  a con- 
strained system),  it  is  readily  appreciated  that  the  cosine  terms  of  real 
and  imaginary  anti  resonances  that  are  far  from  the  forcing  frequency  are 
canceled  by  the  cosine  terms  of  the  natural  frequencies  that  are  also  far 
from  the  forcing  frequency.  The  antiresonant  cosine  terms  near  the 
forcing  frequency  are  very  nearly  one  for  imaginary  anti  resonances. 

Figure  6 shows  a plot  of  the  cosine  terms  for  real  anti  resonances  and 
natural  frequencies  versus  the  nondimensional  forcing  frequencies  and 
damping  coefficients.  It  will  be  noted  that  natural  frequencies  and 
anti  resonances  have  significant  effects  on  the  cosine  terms  only  in  the 
vicinity  of  the  forcing  frequency.  That  is,  the  magnitude  of  the  undamped 
mobility  differs  from  that  of  the  absolute  mobility  only  to  the  extent 
that  the  cosine  terms  differ  from  unity. 

These  cosine  terms  are  cosines  of  v and  6,  not  phase  angles.  The  angles 
v and  6 are  not  physically  realizable  quantities  in  themselves,  as  is  a 
phase  angle,  but  are  the  angles  of  numerator  vectors  of  Equation  77, 

' w2/a2)2  + 92  , relative  to  the  undamped  axis  of  frequency,  and 

similarly,  for  the  denominator.  Although  not  directly  measurable,  they 
are  easily  calculated  from  test  data. 

"R 

Consider  a plot  of  Y.  , which  is  the  real  part  of  3yH/3f  - the  mobility  - 

J A J A 

**R  ** 

and  of  Zjk,  which  is  the  real  part  of  l/(3yj/3fk)  - the  system  impedance. 


i i 


Here  the  system  is  real  impedance,  the  real  part  of  the  reciprocal  of  the 
complex  mobility.  Undamped  anti  resonances  are  obtained  from  the 
impedance  plot  by  the  method  of  Equation  43  for  natural  frequencies. 

The  damping  is  found  by  the  techniques  of  Equation  44.  It  is  extremely 
important  to  note  that  Equation  43  can  be  applied  to  the  impedance  to 
obtain  only  the  real  anti  resonances  and  that  only  the  real  anti  resonances 
are  required  in  this  method. 

Having  found  the  real  anti  resonances,  we  may  form  for  each  real  anti- 
resonance and  forcing  frequency  the  cosine  of  v by  Equation  80.  Similarly, 
the  cosines  of  6 are  formed  from: 
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Figure  6.  Cosines  of  Antiresonant  and  Natural  Frequency  Vector  Angles 


Following  Equation  79,  the  undamped  mobility  at  any  frequency  u)  is 
given  bv 


(83) 


where  the  range  from  i=l  to  N and  i=l  to  M are  generally  in  the  vicinity 
of  the  forcing  frequency.  If  either  M or  N is  odd  and  the  other  even. 


the  proper  sign  of  the  undamped  mobility  may  be  discerned  by  inspection 
of  YR.  Let  cos' 6^  = 

and  substitute  this  for  cos  6-^in  Equation  83.  It  is  seen  then  that 

Equation  83  approaches  Equation  76,  giving  the  modal  acceleration, 
as  a)  approaches  ft.. : 
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At  the  first  real  mobility  peak  for  any  mode,  we  can  calculate 

rtcos  v-p, 
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and  similarly  for  the  p?  real  mobility  peak  for  that  mode.  The  average 
of  the  Y*  for  the  first  real  peak  and  that  for  the  second  will  be  very 
nearly  the  modal  acceleration  of  that  mode.  This  can  be  generalized  into 
one  program  covering  all  modes  in  the  spectrum  of  interest  given  inputs 
of  the  peak  real  acceleration  mobility,  the  frequencies  of  these  peaks, 
the  frequencies  of  the  real  acceleration  impedance  peaks,  the  sign  of 
the  RAC  (E),  and  the  peak  real  acceleration  impedance. 
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LOCAL  SPECTRUM  METHOD 

The  real  mobility  against  frequency  plot  shows  two  peaks  in  the 
vicinity  of  a given  mode.2  Let  the  amplitudes  and  frequencies  of 

these  peaks  be,  respecti  vely,  Yjk  , Yjk  , p]  and  p2  where  pT  > p2< 

Let  the  peak  imaginary  mobility  be  Y*  From  Equations  12  and  13, 
for  the  M-th  mode,  ^ 


where 


yR  _ yR  + y I 

jkp-j  jkp2  'jkM 


'(p2/^)(l  - P f/n^) 

(1  - P-j/^)  + 9m 


S Ajki 
1=1 
i7M 


(puffin  - p2/«m} 

o - p|/ng)2  + 9^ 


(P-I/Qi)2(i  - P?/aj) 

(l  - pf/nf)2  + g? 


(p2/a1)2(i  - p2/^2) 

(1  - p2/ft2)2  + g 2 


9i  (^|/^*j ) 

(1  - fi2/^)2  + g2 


P?/«n  = ]-0+  p2/p|) 


1 


Flannelly  and  Giansante. 
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P2/aM  ’ I <’  + Pf/Pf) 


2 2 
n _ P1  ' P2 

M " 2plp2 

Substitute  Equations  85,  86,  and  87  into  the  terms  not  under  the 
summation  sign  of  Equation  84,  and  drop  the  terms  under  the  summation 
sign  as  they  will  later  be  shown  to  be  negligible.  Define 

y’R  _ yR  . “R 
jkp1  jkp2  AYjk 


then. 


••n  ••  T 

4Vjk  * Ya 


= -A,. 


0 + Pi/PgHi  - pf/Pg)  (p,  + - Pg) 


(1  - p2/p2)2  + 4g 


/ n2  _2 \ **  Pi  9 9 t 

(pi  • f?1  *4-  (p?  - P|) 

Po 


1 + P2  / 1 + 1 \ _ X 

2 ‘ P?  V + P2/Pf  1 + Pf/P2  j 9M 


"'AjkH  p* 


But  i ~ 2 , 2 + ~ ? 2 = Therefore,  using  this  result  and  Equation 

1 + P2/P1  1 + P^/P2 

87  in  Equation  88, 


< * * *J 


(1  + P2/Pt  )' 
M * Ajk"  1 - (p2/p,): 


' ■-  * **  • ’■»'  ^5  1 

assuming  terms  under  the  summation  sign  of  Equation  84  are  negligible. 


(1  - P2/P,) 
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Next,  consider  the  terms  inside  the  simulation  sign  of  Equation  84,  and 
designate  their  sum  as  T.  Let  and  note  that  for  g « g, 

it  can  be  shown  from  Equations  85,  86  and  87  that  p^/  ^ = (1  + gM) 
and  pj;/^  = (1  - gM).  Then: 


T = 


N 

£ A*< 
i^M 


M1  + gH)n  - gjO  + 9M)3 

[l  - e^l  + gM)]2  + g2 


M1  - ^ - M1  - 

[1  - ^.(1  - gM)]2  + g2 


g.g. 

ial 


0 - f^)2  + g2 


(91) 


, O 

It  is  convenient  to  represent  g..  = g^  = g.  For  1 » g » g and 
2 2 

(1  - S)  » g , it  is  seen  that  T « 1 and  the  terms  are  negligible. 
ELECTRONIC  MASS 


Discussion 


This  section  describes  the  theory  of  the  Kaman  Electronic  Mass  (E-MASS) 
and  its  application  to  improved  vibration  testing.  The  E-MASS  is  used 
in  conjunction  with  all  the  information  derivable  from  impedance  tests 
with  speed  and  accuracy  heretofore  unobtainable.  The  E-MASS  method  of 
testing  simplifies  the  calibration  of  impedance  measuring  equipment, 
requires  only  the  coincident  response,  and  eliminates  the  measurement  of 
all  quantities  except  frequency.  E-MASS  is  the  only  method  capable  of 
obtaining  undamped  mobility  at  any  frequency  in  the  spectrum  for  any 
specified  mobility  response.  The  E-MASS  is  expected  to  contribute  sig- 
nificantly to  dramatic  reductions  in  the  developmental  costs  of  new 
helicopters,  airplanes,  and  spacecraft. 

The  E-MASS  electronically  creates  the  dynamic  effects  of  positive  or 
negative  mass  between  the  excitation  and  any  response  point  of  a structure 
in  real  time,  a situation  not  physically  realizable  with  actual  weights. 
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The  E-MASS  device  and  its  associated  software  enable  the  engineer  to 
determine  all  three  parameters  necessary  and  sufficient  to  define  the 
true  mobility  response  of  a structure: 

(1)  The  true  (undamped)  natural  frequencies  used  in  analysis 

(2)  The  structural  damping  mode-by-mode  for  each  mobility 
response 

(3)  The  modal  accelerations  for  each  mode  and  mobility  response 

The  ratios  of  modal  accelerations  determine  the  elements  in  the  modal 
eigenvectors,  and  the  reciprocal  of  a driving  point  modal  acceleration  is 
the  "generalized  mass"  referred  to  that  coordinate. 2 An  accurate 
knowledge  of  the  three  parameters  listed  above  enables  the  engineer  to 
determine,  at  any  point  on  the  structure  he  places  an  accelerometer,  the 
true  complex  response  at  that  point  resulting  from  the  shaking  of  the 
structure  at  any  other  single  point.  He  can  then  perform  system  identifi- 
cation, obtain  static  influence  coefficients  for  arbitrary  restraints, 
determine  maneuver  inertial  loads,  and  accomplish  other  necessary  analyses 
requiring  known  physical  structural  parameters.2*6 

In  actual  practice,  the  engineer  need  only  use  the  E-MASS  system  for 
one  mobility  response,  along  with  conventional  approximations  to  the 
normal  modes  (e.g.  ratios  of  imaginary  mobility  peaks)  if  they  are 
sufficiently  accurate,  to  obtain  all  the  necessary  information.  In 
addition,  E-MASS  can  be  used  at  other  points  to  determine  only  those  modes 
which  are  too  closely  coupled  or  influenced  by  extremely  large  or  ex- 
tremely small  structural  damping  for  accurate  approximation  by  conven- 
tional methods.  Alternatively,  E-MASS  may  be  used  on  all  responses, 
eliminating  conventional  approximations  entirely.  With  the  flexibility 
of  E-MASS,  a variety  of  data  acquisition  methods  can  be  employed  while 
obtaining  all  three  parameters  and  retaining  on-site  data  control  The 
engineer  has  the  opportunity  to  optimize  the  amount  of  information  ob- 
tained per  unit  of  data  recorded  in  a test. 
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Consider  a constant  mobility,  AYjk,  subtracted  from  the  mobility,  Y^. 

••  •• 

The  intersection  of  Y^  and  AY^  describes  new  natural  frequencies  of 

the  system, as  will  be  shown.  Let  be  the  frequency  at  which  AYjk 

intersects  the  undamped  mobility  of  the  system.  These  frequencies  are 
also  the  pseudo-anti  resonances  described  in  Reference  2.  The  new  natural 

frequencies  resulting  from  the  intersection  of  Y..  and  AY.,  are  obtained 
from:  JK 

V - «Jk  <92> 

It  is  mathematically  possible  to  simulate  the  addition  or  subtraction  of 
mass  in  a damped  real  mobility  plot  so  that  the  detrimental  effects  of 
damping  on  the  engineer's  ability  to  analyze  shake  test  data  are  greatly 
lessened.  This  mathematical  idealization  is  physically  meaningless  excepl 
for  the  case  of  a driving  point  mobility.  However,  as  shown  by  the  E-MAS' 
device,  a jk  mass  can  be  electronically  added  or  subtracted  to  a linear 
system  for  producing  new  resonances. 


If  the  impedance  is  changed  by  some  amount  AZ  then. 


V = Z + AZ 

(93) 

or 

Z'  = Y"1  + AZ  = (I  + AZY)Y-1 

(94) 

and  finally 

(ZT1  = Y'  = Y(I  + AZY)"1 

(95) 

If  only  the  jk-th  element  of  the  impedance  matrix  is  changed,  the  new 
jk-th  mobility  is  given  by  Equation  95  as 


(96) 

Anti resonances  are  defined  as  the  frequencies  at  which  the  undamped  mo- 
bility is  zero,  and  natural  frequencies  are  defined  as  the  frequencies  at 
which  the  undamped  mobility  is  infinite.  From  Equation  96, the  addition  of 
a mass  to  the  jk-th  element  changes  the  natural  frequencies  but  does  not 
change  the  anti  resonant  frequencies.  The  new  natural  frequencies  occur 
when  the  denominator  of  Equation  96  goes  to  zero  or 
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Flannel ly  and  Giansante. 


jk 


1 

az: 


jk 


(97) 
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From  Equation  97,  it  is  seen  that  the  new  natural  frequencies  are  exactly 
equal  to  the  frequency  of  the  intersection  of  -l/AZjk  with  the  undamped 

real  mobility.  This  unique  phenomenon  explains  why  the  E-MASS  device 
successfully  determines  the  undamped  modal  accelerations.  The  new  natural 
frequencies  Introduced  by  the  E-MASS  device  are  identical  to  the  pseudo- 
antiresonances described  in  Reference  2. 

Since  the  E-MASS  adds  or  subtracts  mass.  Equation  92  becomes 

4V  - - s (98) 

Combining  Equations  71  and  92  gives 


^jk  Ejk  i?z+]  Ajki  (n2/u>2  _ ^ 


where  AYjk  is  the  E-MASS  setting,  Ejk  is  the  rigid-body  acceleration  co- 
efficient, Ajki  is  the  i-th  modal  acceleration,  is  the  undamped 
resonance,  and  is  the  undamped  pseudo-antiresonance. 

In  matrix  form.  Equation  99  becomes 


AYjk  - 1)]  |Ajkl 


(100) 


The  modal  accelerations  can  be  obtained  independently  of  knowing  the  rigid- 

body  acceleration  coefficient  using  Equation  100  and  two  different 
••  •• 

E-MASS  settings,  AYjk^  and  AYj|<(2)* 


-jk  A*jk(1)j  -i)]|jkl 


(101) 


•jk  ^jk(2)|  [(n2/nZ  . !)]  j jkl 


(102) 


Flannelly  and  Giansante. 
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where  j and  are  the  natural  frequencies  appearing  in  the 

mobility  plot  for  the  corresponding  E-MASS  settings.  Subtracting 
Equation  101  from  Equation  102  gives 


where 


- 4V, 


[Fhi]jAjk( 


(103) 


^ ' <“H)h  - » ' («Hih 


Since  the  rigid  body  acceleration  coefficient  can  be  obtained  from  either 
calculation  or  measurement,  the  accuracy  of  determining  the  modal  accelera- 
tions is  enhanced  by  eliminating  the  rigid-body  acceleration  coefficient 

••  •• 

from  the  calculations  required  by  Equation  100.  If  AY..  and  AY.. 

Jk(l)  Jk(2) 

are  of  opposite  sign,  the  frequency  expressions  for  the  diagonal  terms  of 
Fhi-  add  while  those  for  the  off-diagonal  terms  subtract.  This  mathemati- 
cal feature  serves  to  make  Fhi  more  diagonal,  and  consequently,  better 

conditioned  than  the  matrix  in  Equation  100.  The  E-MASS  settings  and 
the  resulting  shifts  of  the  natural  frequencies  determine  the  modal 
accelerations  from  Equation  103. 


, = N 


(104) 


It  is  not  necessary  to  have  the  same  E-MASS  setting  throughout  the  entire 
frequency  spectrum  since  only  frequency  shifts  are  of  interest.  A change 
in  the  E-MASS  setting,  mode-by-mode,  gives  larger  or  smaller  frequency 
shifts,  allowing  for  more  efficient  and  convenient  operation  of  the 
E-MASS  device. 


57 


The  E-MASS  device  and  Equation  104  are  used  to  determine  for  specified 
response  and  excitation  coordinates  the  one  parameter  not  obtained  by 
present  state-of-the-art  testing  methods. 


The  preceding  development  and  discussion  have  concentrated  on  the  measure- 
ment of  acceleration  mobility  and  unrestrained  systems.  However,  the 
E-MASS  device  is  equally  applicable  to  restrained  systems.  The  only 
difference  between  restrained  and  unrestrained  structures  is  the  inter- 
pretation of  the  rigid-body  acceleration  coefficient  Ejk.  For  restrained 


systems,  Ej^,  is  zero  and  the  aforementioned  methods  remain  unchanged. 

Since  Equation  103  is  independent  of  the  rigid-body  acceleration  co- 
efficient, it  does  not  matter  whether  Ejk  has  some  finite  value  or  is 

zero,  as  in  the  case  of  restrained  structures. 


E-MASS  Calibration  Procedure 


The  most  critical  and  perhaps  the  most  important  phase  of  any  test  pro- 
gram concerns  the  calibration  procedure  necessary  for  relating  electrical 
responses  to  physical ly-meaningful  responses.  The  E-MASS  device  can  be 
used  to  eliminate  most  of  the  work  now  required  for  calibrating  mobility 
response  data  in  terms  of  g/lb,  since  the  E-MASS  output  is  physically 
calibrated  as  a direct  measure  of  mass  ( 1 b/g ) . The  calibration  of  E-MASS 
and  its  applications  in  system  calibration  are  necessary  here. 

The  physical  calibration  of  the  E-MASS  device  is  dependent  upon  a known 
input  condition.  A static  mass  of  known  magnitude  is  mounted  on  an 
impedance  head  of  known  calibration  in  terms  of  lb/volt  and  g/volt  and 
coupled  to  an  electronic  exciter. 

The  outputs  of  the  E-MASS  include  a force  signal  of  the  form 


fout  = f1n  i ^ C05) 

• • 

and  an  acceleration,  y.  The  E-MASS  real  impedance  response 

”R  ^in 

AZ*  = + Am  (106) 

y 

is  obtained  from  the  automatic  mechanical  impedance  analysis  system  and 
recorded  on  the  analog  X-Y  recorder,  yielding  a calibration  line  for  the 
E-MASS  device  in  terms  of  lb/g  versus  the  E-MASS  setting,  AS.  A typical 
plot  of  the  E-MASS  impedance  response  is  illustrated  by  Figure  7. 
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Figure  7.  E-MASS  Impedance  Response 


where 

ZR  = M - Am  = M - K£AS  (107) 

To  facilitate  a general  calibration  constant  for  the  E-MASS  device,  say 
Cemass*  consicler  the  following  development. 

Am(lb/g)  = KeAS  (108) 

where  KE  is  the  slope  of  the  E-MASS  calibration  line  for  known  transducer 
calibrations.  In 

Mlb/vf) 

AmOb/g)  = \C(-q/v")  Am(vf/va)  (109) 

A a 

Kp  and  are  the  force  and  acceleration  transducer  calibration  constants 

respectively,  obtained  from  specified  sensitivity  and  range  settings  on 
the  charge  amplifiers.  Therefore, 


M 

s \ 


and 


where 


Am(vf/va) 


Am(lb/g) 


CEMASS  = 


(110) 

(111) 

(112) 


The  force  and  acceleration  transducers  are  calibrated  through  their 
respective  charge  amplifiers  to  measure  1 lb/volt  and  1 g/volt  for 
specified  sensitivity  and  range  settings.  The  E-MASS  unit  described 
was  calibrated  using  a 2.127-lb  weight.  Figure  8 illustrates  this 
E-MASS  calibration  line,  which  directly  relates  the  amount  of  mass 
electrically  added  or  subtracted  by  the  E-MASS.  A calibration  constant 

°f  CEMASS  = . 00009  7404 ( vf / va ) / AS 

was  obtained  for  this  particular  E-MASS  unit,  which  can  electrically  add 
or  subtract  a maximum  of  50,000  pounds. 

Transducer  Calibration 


Force  and  acceleration  transducers  are  usually  calibrated  in  terms  of 
pi co-coulomb/ lb  or  pcb/g.  The  sensitivity  and  range  settings  on  most 
charge  amplifiers  enable  the  user  to  relate  the  transducer  outputs  as 
lb/volt  or  g/volt.  The  accelerometers  can  be  appropriately  calibrated 
to  match  a "master"  accelerometer  in  terms  of  g/volt.  The  accelerom- 
eters are  matched  by  selecting  a specified  charge-amplifier  range  setting 
and  adjusting  the  charge-amplifier  sensitivity  control  until  the  required 
g/volt  is  obtained.  The  recorded  range  and  sensitivity  settings  serve 
to  calibrate  the  transducers. 

For  example,  suppose  a "master"  accelerometer  is  calibrated  to  measure 
lg/volt.  For  an  acceleration  of  1 g (using  the  "master"  accelerometer) 
and  a specified  charge  amplifier  range  setting,  the  sensitivity  is  ad- 
justed until  the  filtered  accelerometer  signal  measures  .707  volt  rms. 
Similarly,  if  a 2-pound  weight  is  excited  at  1 g,  then  the  sensitivity  is 
adjusted  until  the  filtered  force  signal  measures  1.414  volts  rms.  The 
force  transducer  is  therefore  calibrated  to  measure  1 lb/volt.  The 
transducers  used  to  calibrate  the  E-MASS  device  have  the  following 
characteristics  for  a 2.127-lb  weight: 
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CHARGE  AMPLIFIER 

SETTINGS  AND  OUTPUT 

Transducer 

Range 

Sensitivity 

Output 

Force 

1000 

1.02 

1 lb/ volt 

Accelerometer 

10 

1.058 

1 g/volt 

The  output  of  the  force-charge  amplifier  measures  1.504  volts  rms,  which 
is  equivalent  to  the  2.127-lb  weight.  The  E-MASS  device  was  calibrated 
for  range  settings  of 

Rp  = 1000  = 1 lb/ volt 
RA  = 20  =2  g/volt 

on  the  charge  amplifiers.  For  these  particular  settings,  Kp  is  obtained 
from  Figure  8 and  CEMAS<;  is  known. 

Observe  that  changing  the  charge-amplifier  range  settings  also  changes 
the  physical  slope  of  the  E-MASS  calibration  line.  The  range  settings 
can  be  used  to  relate  Am  to  AS  by 


Am(lb/g)  =(R^J  CEMASSAS 


(113) 


where 


'EMASS  ' 


V «.,  Ke 


(114) 


If  (Ra/M  changes,  then  KF  automatically  compensates  to  give  the  same 
r cal  11 

E-MASS  calibration  constant,  CpMAS<.. 

This  behavior  of  the  impedance  versus  E-MASS  setting  is  illustrated  as 
follows  (Figure  9): 
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Figure  9.  Impedance  Vs  E-MASS  Setting 


CEMASS  = RA/RF 


CEMASS  = 2RA/RF  2AS 


System  Calibration 


The  calibrated  E-MASS  device  serves  as  an  excellent  system  calibration 
method  for  shake  testing  of  restrained  and  unrestrained  structures. 

The  jk-th  acceleration  mobility  defined  by  the  E-MASS  output  is 


f + Amy 


(115) 


If  the  true  force  signal,  f,  is  nulled  for  any  convenient  E-MASS  setting, 
Am,  then  the  calibration  of  the  measured  acceleration  mobility  simply 
becomes 


<vcaris(9/lb) 


(116) 


This  procedure  is  sufficient  to  calibrate  the  entire  mobility  response 
for  known  charge  amplifier  settings  and  transducer  calibrations. 


E-MASS  Data  Reduction 


The  necessary  analytical  expressions  for  determining  the  modal  parameters 
of  a linear,  structurally  damped  system  are  summarized  here. 

Natural  Frequencies: 


(117) 


Structural  Damping: 


Modal  Accelerations: 


(118) 


where  p-j  and  p2  are  the  frequency  peaks  which  occur  in  the  real  damping 

mobility.  Recall  that  and  are  the  undamped  resonances  created 

when  mass  is  electrically  added  or  subtracted  by  the  E-MASS  settings  of 
••  •• 

AYj.  and  AY..  , respectively. 

Jk(2)  jk(1) 


64 


The  application  of  the  E-MASS  in  vibration  testing  is  summarized  as 
f ol 1 ows : 

(1)  To  obtain  the  desired  acceleration  mobility  response  with 
a zero  E-MASS  setting. 

(2)  For  any  convenient  E-MASS  setting,  to  null  the  force  charge 
amplifier  output  and  record  the  pen  position.  The  mobility 
response  is  calibrated  in  terms  of  g/lb. 


Flannel ly  and  Giansante. 


Figure  10  shows  the  test  setup  used  for  the  illustration  of  the  E-MASS. 
Included  in  the  photograph  are  the  beam  specimen,  the  automatic  mechanical 
impedance  analysis  system*  the  E-MASS  device,  and  the  computer  terminal 
used  for  test-site  data  reduction.  Figure  11  isolates  the  E-MASS  device. 
The  beam  specimen,  more  fully  described  in  Reference  2,  was  excited  at 
the  0.5-inch  station,  and  the  response  was  measured  at  the  25-inch 
station.  The  real  damped  transfer  acceleration  mobility  was  recorded 
on  the  analog  X-Y  recorder  and  calibrated  by  the  E-MASS.  The  frequency 
spectrum  of  interest  was  a range  of  10  Hz  to  1200  Hz.  Figures  12  and  13 
illustrate  the  frequency  response  for  the  specimen  with  a zero  setting  on 
the  E-MASS  device.  The  physical  calibration,  frequency  limits,  and  change 
amplifier  settings  are  recorded  on  these  figures.  Figures  14  and  15 
illustrate  the  frequency  responses  for  various  preselected  E-MASS  settings. 
The  necessary  frequency  data  is  read  as  centimeters  from  Figures  14  and 
15  and  converted  to  Hertz  through  computer  software.  The  E-MASS  settings 
are  also  converted  into  g/lb  using  the  E-MASS  calibration  constant  and 
the  charge  amplifier  settings.  Table  1 summarizes  the  natural  frequen- 
cies obtained  for  the  E-MASS  settings  in  g/lb,  obtained  from  the  mobility 
plots.  Equations  117,  118,  and  119  were  used  to  reduce  this  data  to 
obtain  the  natural  frequencies,  modal  accelerations,  and  structural 
damping  of  the  test  specimen.  Table  2 summarizes  the  modal  parameters 
for  the  first  six  modes  of  the  beam.  These  identified  modal  parameters 
were  used  in  conjunction  with  Equations  15  and  16  to  reproduce  the  shake 
test  results.  Figures  16  and  17  present  the  frequency  responses  iden- 
tified by  the  computer  simulation.  Observe  that  the  mobility  amplitude, 
resonances,  and  anti  resonances  are  in  excellent  agreement.  The  E-MASS 
device  can  be  effectively  used  to  obtain  the  necessary  parameters  which 
define  acceleration  mobility  shake  test  data.  Although  this  model  has 
very  little  damping  and  is  perhaps  unrealistic  for  general  structural 
dynamic  applications,  the  E-MASS  device  has  provided  the  engineer  with 
an  excellent  tool  for  analytical  testing.  The  ability  to  use  the  E-MASS 
for  obtaining  appropriate  analytical  models  is  confirmed.  Application  of 
this  analytical  model  to  other  dynamic  problems  has  been  thoroughly  dis- 
cussed in  Reference  2. 

Summary  of  the  E-MASS  Operation 


(3)  To  select  appropriate  intersections  on  the  damped  mobility 
response,  end  to  change  the  E-MASS  setting  until  the  recorder 
matches  the  intersection. 

(4)  To  obtain  the  modified  acceleration  mobilities  for  the  desired 
modes,  repeating  for  each  E-MASS  setting. 

(5)  To  record  the  E-MASS  settings  and  the  corresponding  frequency 
peaks.  To  input  this  data  into  the  software  reduction  system 
and  obtain  the  modal  parameters. 


. 
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Acceleration  Mobility  for  a Response  at  the  25-Inch  Station  and  a Force  a 
0.5-Inch  Station  (Frequency  Range  299.3  - 1224.6  Hertz,  E-MASS  Settinq  is 
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Figure  17.  Acceleration  Mobility  for  a Response  at  the  25-Inch  Station  and  a Force  at  the 
0.5-Inch  Station  (Frequency  Range  299.3  - 1224.6  Hertz,  Response  Obtained  From 
E-MASS  Software) 
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TABLE  1 


NATURAL  FREQUENCIES  OBTAINED  FROM  POSITIVE 
AND  NEGATIVE  E-MASS  SETTINGS 


E-MASS  Setting 


Resonances 


Mode 

As 

In/Sec  /Lb 

_ . (Hz) 

I 

0 

0 

65.4 

II 

0 

0 

187.54 

III 

0 

0 

415.25 

IV 

0 

0 

569.72 

V 

0 

0 

672.52 

VI 

0 

0 

948.95 

I 

115 

68.99 

75.28 

II 

115 

68.99 

245.05 

III 

47 

168.8 

436.83 

IV 

47 

168.8 

546.24 

V 

47 

168.8 

612.44 

VI 

47 

168.8 

894.07 

I 

-115 

-68.99 

50.87 

II 

-115 

-68.99 

158.57 

III 

-140 

-56.67 

400. 75 

IV 

-140 

-56.67 

580.14 

V 

-140 

-56.67 

760.15 

VI 

-47 

-168.8 

1075.9 

m 

I 


r i 
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TABLE  2.  SUMMARY  OF  MODAL  PARAMETERS  FOR  A FORCE  AT 
0.5- INCH  STATION  AND  A RESPONSE  AT  25- INCH 
STATION 


Natural 

Frequencies 


Modal 

Accelerations 


Structural 


DYNAMIC  MODEL  FORCE  DETERMINATION  TESTS 


DYNAMIC  MODEL  TEST  FACILITY 


Description  of  the  Helicopter  Dynamic  Model 


An  experiment  was  undertaken  to  demonstrate  the  ability  of  force 
determination  to  determine  a helicopter's  hub  loads  directly  from 
measured  fuselage  accelerations.  Further,  an  investigation  of  the 
sensitivity  of  force  determination  to  typical  helicopter  gross  weight 
and  center  of  gravity  changes  was  desired.  The  principal  objective  of 
the  test  specimen  was  to  have  a dynamic  model  that  realistically 
simulates  the  characteristics  of  typical  helicopters.  The  dynamic 
characteristics,  in  order  of  importance,  are:  (1)  an  appropriate 
mass  and  stiffness  distribution  - a relatively  stiff  cockpit  and  cargo 
section  with  an  elastically  softer  tail  boom  section;  (2)  a relatively 
low  fundamental  natural  frequency;  and  (3)  representative  structural 
damping.  As  shown  in  Figure  18,  a 5-foot-long  dynamic  model  was 
designed  to  reflect  the  overall  structural  and  mass  characteristics 
of  military  utility  helicopters.  The  model,  constructed  of  riveted 
and  bolted  aluminum,  has  a large  central  fuselage  with  open  sides,  a 
longitudinal  carry- through  structure  of  small  depth  at  the  deck  and 
roof,  and  "beefy"  vertical  members  at  the  forward  and  aft  ends  of  the 
main  cabin.  The  flight  deck  is  essentially  a cantilevered  keel,  and 
the  tail  boom  is  a "bent"  cantilever  beam. 

Steel  weights  (confined  to  a space  of  small  depth)  on  the  roof  structure 
of  the  dynamic  model  simulate  the  masses  of  engines  and  transmission, 
while  steel  and  aluminum  weights  simulate  the  tail  rotor  system  and 
weight  concentrations  on  the  flight  deck,  as  shown  in  Figure  18.  Fuel 
and  payload  were  simulated  by  layers  of  steel  weights  in  the  center  cabin, 
designed  to  be  removable  to  examine  the  effects  of  fuel  burnoff  and  center 
of  gravity  shifts. 

The  first  natural  frequency  of  the  dynamic  model  was  made  to  be  18.1  Hz, 
in  imitation  of  dynamic  data  showing  18  Hz  on  the  Kaman  HU2K  Navy  utility 
helicopter.  Fore-and-aft  adjustments  of  a flight  deck  mass  and  a tail 
boom  mass  allowed  fine  tuning  of  dynamic  model  natural  frequencies.  The 
measured  natural  frequencies  of  the  dynamic  model  were  18.1,  19.6,  40.8, 
65.9,  75.8  and  77.9  Hz.  It  will  be  noted  that  some  of  these  modes  are 
quite  close  together,  reflecting  a situation  not  uncommon  in  actual 
aircraft. 

A unit  weight  and  inertia  breakdown  of  the  model  is  shown  in  Table  3. 

Butt  Line  Zero  is  the  plane  of  symmetry,  the  stations  are  measured  from 
the  nose,  and  the  water  lines  are  measured  from  the  cargo  deck  area. 

These  reference  lines  are  shown  in  Figure  18. 
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PROFILE  VIEW 

Figure  18.  Schematic  of  Dynamic  Test  Model 


TABLE  3.  UNIT  WEIGHTS  AND  MOMENTS  OF  INERTIA  OF  THE  DYNAMIC  MODEL 


Location 


Moment  or  Inertia 


; 

I 


I 


r-- 


■ 


The  total  weight  is  47.78  pounds. 


The  moments  of  inertia  and  products 

1362; 


of  inertia  about  the  center  of  gravity  are  I = 11,280;  I 

xx  yy 

Izz  = 10,391;  Uxy  = -5.77; 


U = -.839;  and  U = 1445  lb-in.  Butt 

At  Jr  £i 

Station  zero  is  the 


plane  zero  is  the  plane  of  symmetry  of  the  model, 
nose,  and  waterline  zero  is  the  cargo  deck. 


Suspension  System 


The  dynamic  model  was  suspended  upside  down,  allowing  a convenient  hub 
shaker  location,  on  bungee  cords  of  a low  spring  rate.  The  shaker  is 
much  larger  relative  to  the  dynamic  model  than  a full-scale  shaker  would 
be  relative  to  an  actual  fuselage.  This  made  it  impossible  to  hang  the 
model  right-side-up  with  a realistic  fore-and-aft  spread  of  bungees. 
Further,  the  upside-down  suspension  simplified  the  work  involved  in  the 
many  shaker  positions  tested,  thereby  reducing  costs. 

A winch  system  for  regulating  the  free  length  of  the  fore-and-aft  bungee 
cords  allowed  the  dynamic  model  to  be  trimmed  in  pitch  during  center  of 
gravity  changes.  It  also  provided  a simple  fine  adjustment  of  the  height 
of  the  dynamic  model,  allowing  the  shaker  armature  to  be  quickly  and 
easily  positioned  between  the  limit  stops  as  the  gross  weight  was  changed. 
The  winch  system  was  used  to  elevate  the  model  out  of  the  way  of  the  work 
table  so  that  calibrations  on  the  impedance-head,  shaker  checks,  and 
other  work  not  involving  the  dynamic  model  itself,  could  be  done. 

The  inverted  "U"  support  of  the  winch  suspension  system  was  attached  to 
the  work  table  and  was  of  adjustable  height  so  that  the  entire  test  set- 
up was  independent  of  any  fixed  structure  and  could  be  moved  for  demon- 
stration purposes  with  little  adjustment  or  dismantling  effort. 

Such  design  features  make  the  test  setup  and  the  dynamic  model  a con- 
venient facility  for  inexpensive  dynamic-model  testing  beyond  the  work 
done  under  this  contract. 

Instrumentation  Provisions  on  the  Dynamic  Model 


As  shown  in  Figure  18,  the  dynamic  model  was  fitted  with  14  accelerom- 
eters, seven  vertically  and  seven  laterally.  Concentrated  weights  on 
the  model  were  drilled  and  tapped  for  mounting  an  impedance  head  later- 
ally or  vertically  at  positions  the  engineer  might  choose.  The  flat 
horizontal  and  vertical  surfaces  of  the  dynamic  model  allow  an  infinite 
number  of  positions  at  which  additional  fixed  accelerometers  can  be 
bonded  or  at  which  a roving  accelerometer  can  be  wax  mounted.  The  use 
of  a roving  accelerometer  in  dynamic  model  testing  must  be  done  with  care 
to  account  for  the  mass  of  the  accelerometer  itself,  as  it  is  not  scaled 
down  from  fwl 1 size. 


The  stud-mounted  fixed  accelerometers  were  mounted  with  special  pre- 
cautions to  avoid  case  distortion,  which  can  produce  physically 
meaningless  signals.  Great  care  was  taken  to  assure  the  surface 
flatness,  the  verticality  of  the  tapped  hole,  and  the  mounting  torque, 
and  that  the  stud  design  met  the  manufacturer's  exacting  specifications. 

The  connection  of  the  impedance  head  from  the  dynamic  model  to  the  shaker 
had  to  be  axially  rigid  but  flexible  in  rotation  and  in  the  lateral 
directions  so  that  restraints  would  not  be  introduced.  The  typical 
necked-down  rod  used  in  testing  larger  specimens  is  not  suitable  for 
scale  model  use,  so  a special  connector  consisting  of  orthogonal  flat 
steel  columns  in  series  was  designed  for  use  on  the  dynamic  model.  The 
plates  of  the  new  sting  were  proportioned  to  carry  axial  loads  much 
higher  than  required  for  shaking  to  preclude  buckling  and  to  have 
virtually  zero  restraint  in  torsion  and  in  the  slopes  and  deflections 
in  the  plane  perpendicular  to  the  shaking  axis. 

The  positions  of  the  accelerometers  on  the  model  are  shown  in  Table  4. 
Odd-numbered  accelerometers  sense  vertical  vibration;  even-numbered 
accelerometers  sense  lateral  vibration. 

With  x,  y and  z measured  from  the  center  of  gravity,  the  RACs  are 
determined  from  Equation  (66). 

For  j and  k odd-numbered  accelerometers, 

Ejk  = 8*0787  + .3323  XjXk  + 1.7349  x 10~4(xjyk  + x^yj) 

+ .03422  y^yk  in. /lb-sec2 

For  j and  k even  numbered  accelerometers, 

Ejk  = 8.0787  + .04357  y.yk  - .0462(y..zk  + y^)  + .3323  zjZ|c 

in. /lb-sec2 


For  accelerometers  in  which  the  indices  sum,  j + k,  to  an  odd  number, 

Ejk  = zj(-1*7349  x 10-4yk  - .3323xk)  + yj(.04621xk  + 2.6889  x 10“5yk) 

in/lb-sec2 
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TABLE  4.  ACCELEROMETER  LOCATIONS 


From  Station 
Zero 
(in.) 


From  Waterline 
Zero 


Accelerometer 


Number 


1 


r 


i 

t ! 
■''.j 


The  damping  characteristics  of  the  dynamic  model  varied  somewhat  from 
those  desired.  However,  the  variations  produced  data  on  the  conserva- 
tive side,  so  that  no  significant  degradation  of  the  experiments  resulted. 

Structural  damping  coefficients  in  the  3-percent  range  were  desired,  but 
in  spite  of  riveted  joints,  the  coefficients  were  in  the  1/2-percent  to 
1-1/2-percent  range.  Consideration  was  given  to  increasing  the  struc- 
tural damping  of  the  dynamic  model  by  using  thin  rubber  laminations  and 
high-hysteresis  coatings,  but  it  was  finally  decided  to  leave  the  damping 
low  and  to  err  on  the  conservative  side.  Low  damping  is  conservative 
because  the  lower  the  damping,  the  steeper  the  mobility  curves  near  res- 
onances, and  hence,  slight  variations  in  the  forcing  frequency  in  these 
spectral  regions  result  in  large  changes  in  mobility.  Therefore,  the 
calibration  error  is  greater  near  resonances.  With  more  damping,  the 
calibration  errors  near  resonances  would  be  smaller,  giving  greater 
accuracy  in  force  determination.  If  the  damping  treatment  of  the  dynamic 
model  were  to  increase  the  damping  to  a level  greater  than  that  which 
would  be  experienced  on  a full-size  aircraft,  the  accuracy  of  force 
determination  shown  in  these  experiments  would  be  unrealistically  good. 

It  was,  therefore,  conservative  to  keep  the  model  damping  on  the  low  side. 

Originally  the  t^il  boom  was  constructed  of  bolted,  aluminum  laminations. 
This  boom  failed! to  give  the  desired  level  and  type  of  damping;  the 
laminations  produced  ankylosis  instead.  Ankylosis  (Reference  12)  is  a 
condition  of  "sticking"  and  "unsticking"  which  sometimes  occurs  in  small- 
scale  models  with  -fetolted  structure  of  the  degrees  of  freedom.  It  was 
first  noticed  in  wind  tunnel  experiments  at  the  National  Physics 
Laboratory  in  England  in  1933  (Reference  12).  In  that  instance,  the 
ankylotic  effect  of  friction  made  an  otherwise  stable  airplane  wind 
tunnel  model  unstable.  This  phenomenon  is  apparently  confined  to  models 
since  there  are  no  known  references  to  its  occurring  in  actual  aircraft. 

Ankylosis  in  this  experiment  was  corrected  by  replacing  the  bolted, 
laminated  tail  boom  with  a tail  boom  of  solid  cross-section.  With 
the  solid  tail  boom,  the  model  was  more  representative  of  a helicopter 
with  a typical  dynamically  well-behaved  semi-monocoque  closed  tail 
section. 

Some  in-house  research  at  Kaman  on  nonmetal  lie  damping,  which  does  not 
bear  on  the  work  in  this  contract,  indicates  that  the  hysteretic  damping 
levels  of  the  dynamic  model  could  be  increased  by  controlled  amounts, 
should  that  be  desired  in  the  future  use  of  the  model.  Its  modular  con- 
struction would  also  allow  the  replacement  of  major  portions  of  the 
structure  with  units  of  composite  material  with  various  damping  and 
dynamic  characteristics. 


matrix  [Y],  The  first  column  of  this  mobility  matrix  is  the  complex 
ratios  of  the  accelerations  at  each  of  the  fourteen  fuselage  accelerom- 
eters to  unit  vertical  forces  at  the  hub.  The  second  column  of  the 
mobility  matrix  is  the  complex  ratios  of  the  accelerations  at  each  of 
the  fourteen  fuselage  accelerometers  to  unit  lateral  forces  at  the  hub. 
This  mobility  matrix  was  obtained  in  two  ways:  (1)  by  direct  shaking  at 
the  hub  with  only  one  shaker  at  a time,  and  (2)  by  shaking  at  the  nose, 
and  not  at  the  hub,  and  using  modal  accelerations.  The  results  of  the 
tests,  comparing  the  two  mobility  matrices,  are  given  below. 

Next,  the  dynamic  model  helicopter  was  "flown"  by  using  vertical  and 
lateral  shakers  at  the  hub  simultaneously.  Different  flight  speeds  and 
maneuvers  were  represented  by  different  magnitudes  and  phasings  of  the 
shaker  forces  at  the  hub.  For  each  flight  condition,  the  amplitudes  and 
relative  phasings  of  the  fourteen  fuselage  accelerometers  were  recorded. 

Having  obtained  the  mobility  matrix  [Y]  from  the  shake  tests  and  the 

"in-flight"  accelerations  of  the  fourteen  accelerometers  {y } from  "flight", 

the  real  and  imaginary  vertical  and  lateral  forces  at  the  hub  were  cal- 

** 

culated  using  the  complex  pseudoinverse  of  the  mobility,  [Y]  : 

{f)R  + itf}1  = [Y]+{y>  (121) 

The  calculated  magnitudes  and  phase  angles  of  these  forces  were  then 
compared  with  the  magnitudes  and  phase  angles  of  the  forces  measured  on 
the  shakers  in  the  "flight".  This  was  done  for  81  test  cases  covering 
54  combinations  of  vertical  and  lateral  forces  over  phase  angles  ranging 
from  0°  to  180°.  Included  in  these  are  cases  representing  calibration 
errors  due  to  fuel  burnoff  effects  on  the  center  of  gravity  and  the  gross 
weight. 

SUMMARY  OF  FORCE  DETERMINATION  TEST  RESULTS 

Table  5 summarizes  the  average  absolute  errors  found  in  81  tests  of 
force  determination.  The  term  "error"  defines  the  percentage  difference 
between  applied  hub  loads  measured  by  a force  transducer  and  the  hub 
loads  predicted  by  force  determination.  In  general,  the  errors  are 
within  the  achievable  limits  of  accuracy  of  the  measurements  of  the 
applied  loads. 
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TABLE  5.  SUMMARY  OF  FORCE  DETERMINATION 


Configuration 

Average  Absolute 
Percent  Error 
in  Vertical 
Force  Amplitude 

Average  Absolute 
Percent  Error 
in  Lateral 
Force  Amplitude 

Average  Absolute 
Degree  Error 
in 

Phase 

Direct  Hub  Shake  - 

Correct  Calibration 

Nominal 

15  Hz  Calibration 

2.49 

16.61 

3.45 

Nominal 

17  Hz  Calibration 

4.08 

7.69 

3.44 

Nominal 

25  Hz  Calibration 

2.36 

4.47 

4.78 

Mass  Change 
25  Hz  Calibration 

6.18 

1.76 

2.76 

CG  Change 

25  Hz  Calibration 

5.59 

2.74 

4.64 

Modal  Acceleration  Test  - Correct  Calibration 


Nominal 

25  Hz  Calibration  5.19  3.69 

Direct  Hub  Shake  - Flight  Configuration  Different 
From  Calibration  Configuration 


Mass  Change 
Nominal 

25  Hz  Calibration 

24.1 

1.97 

CG  Change 
Nominal 

25  Hz  Calibration 

22.4 

3.02 

The  results  indicate  that  force  determination  of  a configuration  is  a 
precise  method  of  measuring  the  magnitudes  and  phasings  of  hub  forces. 

It  is  indicated  by  the  test  results  that  force  determination  calibrate 
of  a nominally  configured  aircraft  with  vertical  and  lateral  hub  excita 
tions  may  have  errors  of  up  to  24  percent  for  flights  with  60-percent 
fuel  burnoff  and  extreme  center  of  gravity  shift,  which  suggests  that  a 
force  determination  calibration  in  a ground  shake  test  should  be  done 
for  several  increments  of  fuel  burnoff  and  longitudinal  shifts  of  center 
of  gravity. 

Caution  should  be  observed  in  extrapolating  from  the  results  obtained 
with  these  model  tests  to  conclusions  about  the  usefulness  of  force 
determination  on  full-scale  aircraft,  in  spite  of  the  excellent  results 
obtained  by  force  determination  in  flights  of  the  DAVI-modified  UH-1H. 
Many  helicopters  have  more  interaction  among  hub  forces  than  could  be 
investigated  in  this  limited  project,  scale  effects  cause  some  degree 
of  unrealism  in  model  testing,  and  the  model  did  not  simulate  the 
dynamics  of  such  things  as  secondary  structures  or  structural  doors. 

DIRECT  SHAKE  FORCE  DETERMINATION  TEST  RESULTS 

25  Hz  Case  - Table  6 

25  Hz  lies  between  the  second  and  third  modes  of  the  helicopter  dynamic 
model,  which  are  at  19.6  Hz  and  40.8  Hz,  and  represents  a fuselage  which 
is  reasonably  well  detuned  from  resonance  with  blade  passage  frequency. 

17  Hz  Case  - Table  7 

17  Hz  is  only  6.1  percent  below  the  18.1 -Hz  natural  frequency  of  the 
helicopter  dynamic  model  and  represents  a fuselage  that  is  not  well 
detuned  from  resonance  with  blade  passage  frequency. 

15  Hz  Case  - Table  8 


The  16.6-percent  error  in  the  lateral  force  at  15  Hertz  was  of  some 
interest.  Why  did  this  condition  show  more  error  than  the  others? 

Also  we  might  ask  "How  accurate  is  accurate?"  The  6-  or  7-percent  error 
obtained  for  the  other  forces  is  certainly  within  measurement  error  and 
requires  no  further  examination.  However,  the  16.6-percent  error  is  a 
significant  change  (not  necessarily  high)  and  was  explained  by  some 
additional  testing.  Recall  that  the  dynamic  calibration  of  the  model 
was  performed  by  independent  force  application  at  the  hub  to  simulate 
actual  operating  conditions.  The  calibration  matrix  was  obtained  from 
the  acceleration  mobility  data  for  each  accelerometer.  When  hub  forces 
were  simultaneously  applied  to  the  model,  new  hub  boundary  conditions 
were  introduced  by  the  shaker  attachment.  The  forces  are  applied  at  the 
simulated  hub  through  "soft"  flexure  units  that  introduce  a resonance 
of  the  system  near  7.0  Hertz.  If  the  hub  calibration  is  again  performed 
with  the  shakers  connected,  then  the  acceleration  mobilities  (and 
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TABLE  6.  CALIBRATION  AT  NOMINAL  GROSS  WEIGHT  AND  CENTER  OF  GRAVITY 


Average 

Error  4.08%  7.693%  3.44 
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consequently  the  calibration)  are  different.  As  the  forcing  frequency 
becomes  much  greater  than  the  lower  resonance  introduced  by  the  shakers, 
the  acceleration  mobilities  converge  to  the  same  measurement.  One  can 
conclude  that  the  percent  error  for  the  15-Hertz  calibration  can  be 
greater  than  for  the  17-Hertz  calibration,  which,  in  turn,  can  be  greater 
than  the  25-Hertz  calibration.  This  is  confirmed  by  observing  that  the 
errors  in  the  lateral  forces  at  15  Hertz,  17  Hertz,  and  25  Hertz  are 
16.61  percent,  7.69  percent,  and  4.47  percent,  respectively. 

A comparison  of  acceleration  mobilities  with  and  without  the  shakers 
attached  to  the  model  showed  a greater  difference  for  the  lateral  force 
than  for  the  vertical  force.  An  accelerometer  on  the  tail  section  was 
selected  for  comparison  purposes.  This  difference  in  mobility  measure- 
ments explains  why,  at  17  Hertz,  the  vertical  force  calculations  are 
biased  on  the  low  side  while  the  lateral  force  calculations  are  biased 
on  the  high  side.  The  calibrations  for  the  vertical  forces  are  higher 
with  the  shakers  attached,  and  the  calibrations  for  the  lateral  forces 
are  lower.  At  25  Hertz,  this  bias  of  the  calculated  forces  is  not 
nearly  as  noticeable  since  the  effects  of  the  lower  resonance  are  no 
longer  significant. 

The  influence  of  the  shaker  attachments  on  the  method  of  force  determina- 
tion is  primarily  of  academic  interest.  Under  actual  operating  con- 
ditions, this  boundary  condition  effect  will  not  be  present.  The 
agreement  between  the  known  applied  forces  and  the  forces  obtained  via 
force  determination  is  excellent,  and  consequently,  a recalibration  of 
the  model  to  account  for  these  artificial  boundary  conditions  does  not 
appear  beneficial.  It  is  immaterial  whether  the  error  is  7 percent  or 
4 percent  since  these  differences  are  well  within  the  range  of  measure- 
ment error.  Of  course,  if  this  boundary  condition  anomaly  had  been 
worse,  additional  testing  would  have  been  recommended. 

In  order  to  minimize  the  effects  of  the  shaker  attachments,  the  25-Hertz 
condition  was  selected  for  purposes  of  comparison  between  the  nominal 
configuration  and  the  mass  change  and  center- of- gravity  change  configura- 
tions. At  25  Hertz,  the  acceleration  mobilities  of  the  suspended  model 
and  the  acceleration  mobilities  with  the  shakers  attached  appear  to  be 
identical. 

Weight  Change  Case  - Tables  9 and  10 

Any  aircraft  used  in  force  determination  would  be  dynamically  calibrated 
by  ground  shake  tests  for  the  configuration  to  be  flown.  However,  fuel 
burnoff  in  flight  causes  mass  and  vertical  center- of- gravi ty  changes  that 
affect  the  dynamics  of  the  ship.  To  help  evaluate  how  much  fuel  burn- 
off can  be  tolerated  for  a given  calibration,  fuel  burnoff  was  simulated 
by  removing  3.75  pounds  of  steel  weights  from  the  lower  fuselage 
near  the  center-of- gravity  station.  This  represents  a 7.8-percent  change 
in  the  fuselage  weight  (without  rotors).  With  10  percent  of  a helicopter's 
total  gross  weight  being  the  weight  of  the  rotor  and  12  percent  of  a 
typical  helicopter's  gross  weight  being  the  weight  of  the  fuel,  this 
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weight  change  in  the  dynamic  model  represents  about  a 60-percent  fuel 
burnoff. 

Table  10  shows  the  results  obtained  with  the  calibration  corrected  for 
the  change  in  weight  due  to  fuel  burnoff. 

Center- of- Gravity  Change  Case  - Tables  11  and  12 

To  give  an  indication  of  how  much  center-of- gravity  shift  should  be 
allowed  without  calibrating  for  the  shift  in  the  shake  test,  the  center 
of  gravity  of  the  dynamic  model  was  shifted  from  the  nominal  21.52-inch 
station  to  the  20.83-inch  station  for  the  same  gross  weight.  This  3.25- 
percent  change  in  the  center-of-gravity  station  represented  a shift  from 
the  nominal  position  to  one  of  the  extreme  center-of-gravity  positions 
in  a typical  helicopter  and  is  more  severe  than  one  would  actually  allow 
without  accounting  for  it  in  the  shake  test  calibration.  When  the 
corrected  calibration  for  the  center-of-gravity  shift  was  used,  the 
results  of  Table  12  were  obtained. 

MODAL  ACCELERATION  FORCE  DETERMINATION  TEST  RESULTS  - NOSE  SHAKING 
Discussion 

The  purpose  of  these  sets  of  tests  was  to  demonstrate  that  force  deter- 
mination can  be  successfully  accomplished  without  shaking  at  the  conven- 
tional point  of  excitation,  i.e,  the  hub.  With  modal  acceleration 
testing,  the  engineer  can  calibrate  the  fuselage  for  any  number  of  ex- 
ternal forces  by  shaking  at  only  one  point,  although  it  is  advantageous 
to  choose  two  or  three  shaking  points  to  check  the  data's  reliability, 
improve  accuracy,  and  simplify  the  data  reduction. 

The  shaking  point  or  points  for  modal  acceleration  testing  is  a choice 
of  convenience  for  the  engineer.  In  the  interests  of  maximum  accuracy, 
the  shaking  point  or  points  should  have  high  driving  point  mobilities  for 
the  modes  of  consequence,  a condition  which  also  allows  the  use  of 
relatively  small  shakers.  Such  points  are  commonly  near  the  nose  or 
tail  of  a helicopter. 

Procedure 

The  nose  of  the  helicopter  dynamic  model  was  chosen  as  the  shaking  point, 
although  testing  showed  that  the  tail  would  have  been  almost  equally 
suitable.  The  helicopter  dynamic  model  was  shaken  both  laterally  and 
vertically  to  take  advantage  of  the  high  degree  of  decoupling  resulting 
from  symmetry.  The  first  and  third  modes  showed  up  very  strongly  in  the 
vertical  driving  point  mobility  plot  at  the  nose,  while  the  second  mode 
was  negligible.  The  second,  fourth,  fifth,  and  sixth  modes  had  strong 
driving  point  responses  at  the  nose  in  lateral  shake,  while  the  first  and 
third  modes  were  negligible. 
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The  modal  accelerations  for  all  14  accelerometers,  including  the  two  at 
the  hub,  were  determined  using  E-MASS  for  the  vertical  and  lateral  shake. 
The  ratios  of  each  to  its  driving  point  were  compared  to  the  other  as  a 
check  on  the  validity  of  the  data.  For  each  of  the  14  accelerometers, 
the  modal  accelerations  and  RACs  (rigid  body  acceleration  coefficients 
known  from  weight  data)  were  fed  into  a Hewlett-Packard  2000F  time-share 
computer  using  a telephone  line.  The  mobility  curves  calculated  by  the 
computer  using  the  modal  accelerations  and  RACs  were  plotted  on  the  time- 
share  plotting  terminal  and  compared  to  the  actual  test  plot  for  each 
accelerometer  as  a check  on  the  testing.  The  low-frequency  asymptote  of 
the  measured  mobility  curve  on  the  test  x-y  recorder  was  compared  to  the 
RAC,  calculated  from  weights  data,  as  still  another  check.  Agreement  was 
excellent  for  each  case.  These  were  done  for  both  the  vertical  shake  and 
the  lateral  shake  during  the  course  of  the  testing.  In  every  mobility 
plot,  the  hysteretic  damping  coefficient,  g,  and  the  undamped  natural 
frequency  were  determined  for  each  mode. 


The  importance  of  using  the  on-site  data  checks  described  above  was 
demonstrated  early  in  these  tests  by  the  detection  through  these  methods 
of  seriously  distorted  signals  from  two  of  the  14  accelerometers.  With 
conventional  methods  of  testing,  it  is  impossible  to  detect  such  false 
signals.  The  problem  with  the  two  accelerometers  was  corrected, 
and  the  tests  on  the  helicopter  dynamic  model  were  continued. 

Orthonormal  Modes  and  Mobility  Matrix 

The  orthonormal  modal  vector  for  each  mode  was  formed  using  either  the 
modal  accelerations  from  the  nose  vertical  shake  or  those  from  the  nose 
lateral  shake,  depending  on  which  had  the  larger  driving  point  modal 
acceleration.  In  the  first  mode,  for  example,  the  vertical  nose  shake 
had  a driving  point  modal  acceleration  of  11.69  in. /lb-sec2  while  the 
lateral  nose  shake  had  a driving  point  modal  acceleration  of  only  0.17 
in. /(lb-sec2),  so  the  vertical  shake  modal  accelerations  were  used.  The 
nose  vertical  station  was  designated  "1".  The  first  mode  (I)  orthonormal 
vector  for  any  accelerometer  j was  found  from  Equation  20: 


= iiUL 

VAn  i 


in.1/2/(lb1/2-sec) 


As  another  example,  the  lateral  driving  point  modal  acceleration  for 
the  second  mode  was  13.05  in. /(lb-sec2)  while  the  vertical  was  .227 
In. /(lb-sec2),  and  of  course,  the  modal  accelerations  for  lateral  nose 
shake  were  used  In  forming  the  orthonormal  mode  vector  for  the  second 
mode: 
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for  the  j-th  accelerometer  with  nose  lateral  being  designated  as  2. 

The  orthonormal  mode  vectors  for  the  first  six  modes  in  the  helicopter 
dynamic  model  are  shown  in  Table  13. 

Of  the  14  accelerometers,  two  were  positioned  for  hub  excitations,  giving 
a 12x2  matrix  of  mobilities  if  the  hub  accelerometers  were  removed  for 
flight.  In  an  actual  helicopter  one  would  use  as  many  as,  but  not  more 
than,  six  accelerometers  at  the  rotor  hub  during  shake  tests.  The  hub 
accelerometers  would  be  removed  before  flight. 

After  calibration  in  the  above  manner,  the  helicopter  dynamic  model  was 
"flown"  as  described  earlier  with  simultaneous  vertical  and  lateral  hub 
forces  of  varying  magnitudes  and  phasings,  and  the  fuselage  accelerations 
were  recorded.  These  "flight"  accelerations,  premultiplied  by  the  pseudo- 
inverse of  the  mobility  matrix  as  shown  above,  give  the  force  determina- 
tions of  the  magnitudes  and  phases  of  the  forces  applied  to  the  hub  in 
"flight".  The  force  determination  results  are  compared  to  the  actual 
magnitudes  and  phasings  of  the  hub  forces  applied  in  the  "flights"  in 
Table  14. 


TABLE  13.  ORTHONORMAL  MODES 

(in./lb-sec2)t/2 

Mode  Designation  I II  III  IV  V VI 

Natural 

Frequency 

(Hz) 18.1  19.6  46.8  65.9  75.8 77.9 

Accelerometer 


1 

3.419 

.562 

-5.840 

-.199 

-.116 

-.068 

2 

-.537 

3.613 

-1.873 

4.254 

1.683 

2.833 

3 

1.174 

.132 

.429 

0.000 

.088 

-.061 

4 

-.454 

1.574 

.043 

-.320 

-.055 

-.359 

5 

-1.696 

-.329 

.799 

.038 

-.832 

.538 

6 

.171 

-1.837 

-.126 

.081 

-.136 

-.287 

7 

-1.956 

-.045 

-1.840 

-.155 

1.180 

-.694 

8 

.626 

-4.747 

-.327 

2.941 

.027 

.309 

9 

.088 

-.081 

-1.415 

-.039 

3.960 

-3.009 

10 

.058 

-.379 

-.265 

5.075 

-.772 

-5.160 

11 

2.291 

.247 

-.441 

.025 

3.422 

-2.342 

12 

-.188 

1.825 

-.154 

6.327 

-1.455 

-6.888 

13 

6.638 

.883 

1.463 

-1.105 

0.000 

2.109 

14 

-.520 

4.478 

.307 

-3.802 

2.779 

1.936 
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TABLE  14.  MODAL  ACCELERATION  CALIBRATION  AT  NOMINAL  GROSS  WEIGHT  AND  CENTER  OF  GRAVITY 
FLIGHT  AT  NOMINAL  GROSS  WEIGHT  AND  CENTER  OF  GRAVITY 
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Average  5.19%  3.69%  4.40 


The  two  columns  of  the  matrix  of  mobilities  at  frequency  u>  relating  the 
accelerometers  to  the  hub  in  these  dynamic  model  tests  are  given  by: 


II 


for  hub  excitations  designated  5 and  6. 


1 


' 


s 


FL IGHT  APPLICATION  OF  FORCE  DETERMINATION  TO  THE  UH-1 H 


Force  determination  was  used  as  an  adjunct  to  another  project,  the  flight 
evaluation  of  Kaman's  Dynamic  Anti  resonant  Vibration  Isolator  (DAVI). 

This  project  was  also  performed  under  contract  for  the  Eustis  Directorate, 

USAAMRDL.  This  particular  application  of  force  determination  was  nec- 
essarily rather  crude,  but  the  results  were  very  revealing. 

Prior  to  the  installation  of  the  DAVI  rotor  isolation  system,  the  UH-1H 
test  vehicle  was  flown  for  baseline  vibration  data,  and  the  magnitudes 
and  phases  of  accelerations  in  the  fuselage  were  recorded.  The  standard 
UH-1H  was  also  shake  tested  to  obtain  baseline  mobility  data.  Six  ver- 
tical and  six  fore-and-aft  accelerometers  were  used  in  the  baseline  shake 
test  and  flight  test  to  determine  the  two-per-rev  vertical  and  longi- 
tudinal forces  at  the  hub  using  the  force  determination  method.  Pitching 
and  rolling  moments  were  considered  to  be  negligible  because  the  UH-1  has 
a teetering  rotor. 

I 1 

When  the  DAVI  rotor  isolation  system  was  installed,  the  ship  was  shaken 
to  determine  the  mobilities  using  the  direct  shake  technique.  The  hub 
forces  found  through  force  determination  for  the  standard  UH-1H  were 
then  applied  to  the  shake  test  mobilities  on  the  DAVI-modified  ship  to 
predict  the  new  flight  vibrations  with  the  DAVI.  The  comparison  of  these 
force  determination  predictions  to  actual  flight  results  of  the  DAVI  UH-1 
is  shown  in  Figure  19. 

The  force  determination  predictions  correlate  very  well  with  flight  data 
for  the  pilot's  seat  and  the  tail  at  all  flight  speeds  from  20  to  120 
knots.  The  correlations  at  the  pilot's  seat  and  the  nose  are  excellent 
up  to  80  knots.  Beyond  80  knots,  the  force  determination  prediction  and 
the  flight  measurements  diverge.  Force  determination  theory  would  suggest 
that  this  divergence  Is  caused  by  a two-per-rev  vertical  vibratory  force 
acting  on  the  fuselage  In  flight.  The  copilot's  seat  and  nose  curves 
given  by  force  determination  prediction  and  by  flight  measurement  indicate 
that  this  unknown,  two-per-rev  fuselage  force  becomes  effective  at  about 
80  knots  and  Is  a function  of  the  forward  speed  of  the  ship.  This 
suggests  an  aerodynamic  excitation  other  than  downwash. 

These  predictions  from  force  determination  led  to  the  suspicion  that  the 
horizontal  stabilizer,  mounted  on  the  tall  boom  and  actuated  by  motions 
of  the  swashplate,  was  the  excitation  source.  There  Is  a two-per-rev 
relative  motion  on  the  horizontal  stabilizer  pitch  control  in  both  the 
DAVI-equipped  aircraft  and  the  conventional  aircraft. 
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Figure  19.  Two-Per-Rev  Vertical  Vibration  on  the  UH-1H 


The  horizontal  stabilizer  was  instrumented  to  check  for  pitch  changes, 
and  test  flights  showed  angle  of  attack  changes  on  the  horizontal 
stabilizer  at  the  two-per-rev  frequency  in  flight,  producing  significant 
two-per-rev  vertical  vibratory  forces  on  the  tail  boom  at  speeds  above 
80  knots.  The  effect  was  found  in  flight  tests  of  both  the  DAVI-equipped 
and  the  conventional  UH-1H. 

The  shake  tests  in  the  DAVI  project  were  done  before  the  force  determina- 
tion method  was  developed.  Had  force  determination  methods  been  avail- 
able at  that  time,  the  unknown  force  could  have  been  revealed  in  the  first 
baseline  flights  by  the  tail  boom  accelerometers. 

This  makeshift  use  of  force  determination  on  an  incidental  and  rather 
limited  basis  not  only  contributed  to  the  success  of  the  DAVI  flight 
test  program  but  demonstrated  the  accuracy  of  force  determination  and 
its  usefulness  in  identifying  unsuspected  sources  of  vibratory  loadings. 
The  broader  and  more  important  implications  of  this  experience  on  future 
dynamic  flight  and  ground  testing  of  helicopters  of  all  types  are  quite 
obvious. 


CONCLUSIONS 


1.  Force  determination  technology  is  sufficiently  established  to  allow 
it  to  be  fully  demonstrated  on  a full-size  helicopter  to  determine 
in-flight  vibratory  main  rotor  loads,  tail  rotor  loads,  aerodynamic 
surface  loads  and  loads  from  external  stores  using  only  fuselage 
accelerometer  data. 

2.  Using  force  determination  flight  test  results,  accelerated  relia- 
bility testing  on  helicopter  airframe  and  equipment  in  ground 
vibration  tests  can  be  done. 

3.  Force  determination  can  be  used  to  correlate  the  rotor  loads  predicted 
by  rotor  analysis  computer  programs  with  the  rotor  loads  and  ship  vi- 
brations actually  measured  in  flight. 

4.  The  methods  of  force  determination  can  be  employed  in  conjunction 
with  rotor  analysis  computer  programs  to  predict  the  vibration 
effects  of  a proposed  new  rotor  on  the  fuselage  and  fuselage-mounted 
equipment,  including  engines  and  gearboxes.  This  may  be  done  using 
flight  test  data  on  an  existing  aircraft  or  analytical  mobility  data 
on  a proposed  new  aircraft  fuselage. 

5.  In  a ground  vibration  survey,  it  is  not  necessary  to  shake  a heli- 
copter at  the  hub  or  to  apply  moments  at  the  hub.  The  data  required 
can  be  acquired  by  shaking  at  one  or  more  convenient  points  using 
modal  acceleration  testing. 

6.  Modal  acceleration  testing  can  be  done  with  relatively  small  shakers 
and  negligible  tare  positioned  at  high  mobility  points,  commonly 
the  nose  and  tail.  This  yields  as  much  information  as  would  shaking 
at  every  accelerometer  location. 

7.  Modal  acceleration  testing  yields  the  classical  analytical  normal 
modes  and  undamped  natural  frequencies  needed  for  practical  com- 
parison to  theory,  as  opposed  to  mere  mobility  ratios  ordinarily 
found  from  shake  tests. 
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LIST  OF  SYMBOLS 


i-th  anti  resonant  frequency 

i-th  anti  resonant  frequency  for  response  at  j due  to 
force  at  x 


modal  acceleration  of  the  i-th  mode  for  response  at  j due 
to  force  at  k - in./secvlb 


E-MASS  calibration  constant  defined  in  Page  60 
jk-th  rigid-body  acceleration  coefficient  - in. /lb-sec 


EMASS 


force  - lb 


structural  damping  coefficient;  also  acceleration  of 
gravity  - in./sec2 


structural  damping  coefficient  of  the  i-th  mode 


imaginary  operator 


moment  of  inertia  - lb-in. -sec' 


stiffness  - lb/in 


accelerometer  transducer  calibration  constant 


slope  of  E-MASS  calibration  curve  - lb/g/As 
force  transducer  constant  - lb/ volt 


generalized  stiffness  of  the  i-th  mode  - lb/in 


mass 


mass 


moment  around  x axis  - lb-in 


generalized  mass  of  the  i-th  mode  - lb-sec  /in 


number  of  modes 


frequency  peaks  defined  in  Page  36 


LIST  OF  SYMBOLS  (Continued) 


\ ; 


r; 

' I 

h 


r 


J 

<1 


i *. 


AS 

U 


z 

z 


♦jl 


•• 

e 

*jl 

U) 


accelerometer  charge  amplifier  range  setting 
force  charge  amplifier  range  setting 
E-MASS  potentiometer  setting 

2 

product  of  inertia  - lb-in. -sec 

output  voltage  of  accelerometer  charge  amplifier  - volts 
output  voltage  of  force  charge  amplifier  - volts 
displacement  - in. 

2 

vibratory  acceleration  - in. /sec 

acceleration  mobility  - in. /lb-sec2 

2 

acceleration  mobility  - in. /lb-sec 

2 

E-MASS  acceleration  mobility  setting  - in. /lb-sec 

2 

displacement  impedance  - lb-sec  /in. 

2 

acceleration  impedance  - lb-sec  /in. 
impedance  change  - lb/in. 
j-th  element  in  the  i-th  normal  mode 
normal  mode 

o i/p 

orthonormal  mode  (in. /lb-sec  ) 

2 

rotational  acceleration  - rad/sec 
j-th  orthonormal  mode  element  for  i-th  mode 
forcing  frequency  - rad/sec,  Hz 
i-th  natural  frequency  - rad/sec,  Hz 
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BRACKETS 

LIST  OF  SYMBOLS  (Continued) 

[ ] 

matrix 

M 

diagonal  matrix 

[i] 

identity  matrix 

{ } 

column  or  row  matrix 

SUPERSCRIPTS 

R 

real 

I 

imaginary 

-1 

Inverse 

T 

transpose 

-T 

Inverse  transpose 

♦ 

pseudoinverse 

* 

a modal  parameter 

SUBSCRIPTS 

1 

modal  Index 

j 

accelerometer  Index 

k 

accelerometer  Index 

X 

accelerometer  Index  forcing  station 

q 

generalized  index 

i, ii  & in 

used  as  subscripts  to  modal  accelerations 

LIST  OF  SYMBOLS  (Continued) 


OTHER  INDICES 

N number  of  degrees  of  freedom 

0 number  of  generalized  coordinates 

HxF  capital  letters  under  matrices  indicate  the  number  of 

rows  and  columns,  respectively 

dots  over  a quantity  indicate  differentiation  with 
respect  to  time  or  a parameter  related  to  differentiation 
with  respect  to  time 

53  summation  symbol  - indicates  the  sum  of  the  indexed  terms 

TT  product  symbol  - indicates  multiplication  of  the  indexed  terms 


